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ABSTRACT 


nils  report  presents  the  results  of  an  analysis  of  the  turbulence  data  derived  from 
the  high  altitude  clear  air  turbulence  (HICAT)  flight  investigation  and  from  the  perti¬ 
nent  meteorological  and  geophysical  data.  The  curves  representing  the  power  spectra 
have  various  slopes  and  shape.  These  curves  are  grouped  into  categories.  One 
group  of  categories  are  predominantly  cases  of  flights  over  mountains.  A  second 
group  of  categories  are  mostly  from  flights  over  flatland.  Theoretical,  statistical, 
and  practical  forecasting  analyses  in  addition  to  flight  records  all  indicate  that  large 
horizontal  temperature  gradients  and  large  changes  in  vertical  temperature  gradients 
are  closely  associated  with  HICAT  and  conversely  small  horizontal  temperature  gra¬ 
dients  and  nearly  constant  vertical  temperature  gradients  are  closely  associated 
with  smooth  air. 

Evidence  is  presented  that  the  relative  frequency  of  occurrence  of  high  altitude  clear 
air  turbulence  by  geographic  area  may  be  deduced  from  an  analysis  of  the  standard 
deviations  of  the  24  hour  temperature  changes  at  the  HICAT  flight  altitude. 

Turbulence  was  observed  during  6.  9%  of  the  flight  miles  over  mountains  and  5. 5% 
of  the  flight  miles  over  water.  The  RMS  of  the  vertical  gust  velocity  during  HICAT 
was  0. 86  over  water,  1. 02  over  flatland  and  1. 89  over  mountains. 

(This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
governments  or  foreign  nationals  may  be  made  only  with  prior  approval  of  the  Air 
Force  Flight  Dynamics  Laboratory  (FDTR),  Wright-Patterson  AFB,  Ohio  45433.) 
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SECTION  I 


INTRODUCTION 


This  report  consists  of  an  extensive  analysis  of  the  meteorological  and  geophysical 
conditions  associated  with  high  altitude  clear  air  turbulence.  The  basic  turbulence 
data  were  obtained  from  Project  HICAT  (Investigation  of  High  Altitude  Critical 
Atmospheric  Turbulence)  of  the  Flight  Dynamics  Laboratory,  Air  Force  Systems 
Command,  Wrlgbt-Patterson  Air  Force  Base,  Ohio. 

The  discussion  begins  with  a  description  of  the  basic  turbulence  data  and  of  the 
method  of  selecting  samples  for  analysis.  Then  the  meteorological  and  geophysical 
data  that  were  used  in  the  study  are  described.  Sections  on  case  history  studies, 
statistical  studies,  theoretical  analysis  and  climatological  studies  follow. 

The  purpose  of  fills  report  has  been  to  present  data,  models,  and  summaries  rela¬ 
ting  to  high  altitude  clear  air  turbulence  in  a  form  that  will  provide  a  foundation 
for  the  development  of  Improved  methods  for  forecasting  turbulence,  ejaculating 
gust  loads  and  vehicle  responses,  designing  aerospace  vehicles  to  withstand  gusts, 
and  analyzing  fatigue. 

CHRONOLOGY 

The  HICAT  program,  initiated  in  1962,  was  divided  into  two  distinct  cut  comple¬ 
mentary  phases.  One  phase  consisted  of  the  design,  installation,  and  flight  of 
instruments  that  would  yield  data  relating  to  the  motions  of  file  atmosphere.  The 
first  flight  of  the  WU-2  aircraft  under  this  program  was  made  on  20  February  1964, 
and  by  15  July  1964,  18  HICAT  search  flights  were  completed.  Complete  docu¬ 
mentation  of  this  portion  of  the  HICAT  program  was  presented  by  Crooks  (1). 

New  and  more  sophisticated  instruments  were  then  installed  in  the  WU-2  and  HICAT 
search  flights  were  resumed  on  8  October  1965.  By  1  March  1967,  a  total  of  114 
more  flights  were  made.  The  details  of  file  instrument  development,  data  reduction, 
and  flight  program  for  the  period  ending  in  February  1967  were  presented  by  Crooks, 
Hoblit,  and  Prophet  (2).  A  report  covering  the  flights  made  from  March  1967 
thropjh  February  1968  is  in  preparation. 

The  second  phase  of  the  HICAT  program  was  initiated  in  July  1962.  Work  accomp¬ 
lished  vp  to  April  1965  was  presented  by  Hildreth,  et  al.  (3)  and  Hildreth,  Court, 
and  Abrahms  (4).  The  present  report  covers  the  work  accomplished  up  to  May  1968. 

The  HICAT  program  is  now  part  of  a  larger  ALLCAT  program  on  "Critical  Atmo¬ 
spheric  Turbulence"  which  was  initiated  by  the  Air  Force  Flight  Dynamics  Labor¬ 
atory  in  June  1966.  See  Anonymous  (5).  The. objective  of  the  ALLCAT  Program 
is  to  provide  fine-scale  true  gust  velocity  measurements  of  continuous  wavelength 
turbulence  for  all  altitudes,  seasons,  and  typical  geographic  areas.  These  turbu¬ 
lence  data,  correlated  with  meteorological  eondttons  and  synthesized  with  other 


1 


data,  will  be  used  to  establish  valid  turbulence  criteria  for  the  design  of  advanced 
vehicle  structures  and  flight  control  systems,  and  for  human  factors  investigations. 
The  turbulence  data  can  be  utilized  in  evaluating  structural  modifl  cations  and  life 
determinations  of  existing  flight  vehicles.  These  data  will  assist  meteorological 
agencies  in  verifying,  changing,  and  establishing  techniques  for  fore'jasting  critical 
atmospheric  turbulence.  The  turbulence  data  will  also  be  useful  to  other  technical 
disciplines,  such  ss  physiological  and  psychological  investigations  of  crew  members 
crew  accommodations*  guidance,  flight  instrumentation,  jet  engine  design,  and 
mission  pluming  in  operational  organizations. 


SECTION  n 


OTHER  STUDIES 


Many  reports  and  articles  have  appeared  in  the  last  few  years  on  clear  air  turbulence 
and  its  effects  upon  aircraft.  C!  more  than  300  entiles  in  the  Bibliography  (Section 
V),  many  are  efforts  to  define  the  meteorological  and  geophysical  conditions  associ¬ 
ated  with  dear  air  turbulence  at  altitudes  of  20, 000  to  40, 000  ft.  Observations  of 
clear  air  turbulence  upon  which  these  studies  are  based  vary  from  subjective  reports 
of  airline  and  military  pilots  to  specific  and  detailed  measurements  of  the  character¬ 
istics  of  the  turbulence.  likewise,  the  discussions  vary  from  purely  theoretical 
treatments,  to  case  histories  and  to  extensive  statistical  studies. 

With  few  exceptions,  the  meteorological  data  are  the  temperatures,  pressures,  wind 
velocities,  and  wind  directions  obtained  routinely  by  national  weather  services. 
Horizontal  and  vertical  temperature  gradients  and  wind  shears  may  be  derived 
firom  these  basic  measurements  by  an  analysis  system  that  is  nearly  always  partially 
subjective.  Other  functions  of  the  temperature  and  wind  fields  such  as  Richardson 
number,  deformation,  divergence,  etc.  have  also  been  used  in  statistical  analyses 
and  case  history  studies  of  relationships  between  meteorological  conditions  and  the 
reports  of  turbulence.  In  general,  derivation  of  such  quantities  from  the  basic 
temperature  and  wind  data  requires  a  significant  amount  of  subjective  judgment 
which,  along  with  the  errors  in  the  original  data,  introduces  an  unknown  uncertainty 
in  the  final  results. 

The  basic  meteorological  data  and  the  quantities  derived  from  these  data  that  have 
been  used  in  statistical  studies  of  clear  air  turbulence  are  summarized  in  Table  1. 
These  studies  have  served  as  an  important  guide  to  the  work  described  in  this  report. 


TABLE  1  LITERATURE  SURVEY  OF  RELATIONSHIP  BETWEEN  CLEAR  AIR  TURBULENCE 
AND  METEOROLOGICAL  St  GEOPHYSICAL  PHENOMENA 


Kadlao,  P.W.  (1963)  1960-63  EaatU.S.  Commarolal  80-40.000'  414  fllfhU  4  olrrui  pattam  Inodalai  polar  li  aubtrop- 
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SECTION  m 


TURBULENCE  MEASUREMENTS  AND  THEORIES 


The  case  histories  and  the  statistical  studies  presented  in  this  report  are  based 
directly  upon  the  WU-2  HICAT  flight  program.  From  February  1965  to  March  1968, 
a  total  of  47  hours  of  clear  air  turbulence  were  recorded  in  flights  from  bases  in 
California,  Massachusetts,  Alaska,  Hawaii,  Puerto  Rico,  New  Zealand,  Australia, 
England,  Florida,  Panama,  Maine,  and  Louisiana.  The  flights  and  measurements 
derived  from  the  flights  from  February  1965  to  March  1967  are  discussed  in  detail 
by  Crooks,  Hoblit  and  Prophet  (2)  and  this  material  was  available  for  use  in  the 
present  study.  Portions  of  the  turbulence  data  from  the  flights  made  after  March 
1967  were  also  available  for  use  in  the  present  report. 

The  identification  of  turbulent  cases  and  the  derivation  of  the  relevant  in-flight  data 
was  directed  by  Crooks  and  Hoblit.  Because  the  selection  of  turbulence  samples  is 
of  fundamental  importance  for  the  present  study,  the  description  of  the  method  of 
selecting  a  turbulence  sample  will  be  summarized  here. 

The  primary  airborne  data  acquisition  system  is  a  digital  pulse  code  modulation 
system.  The  system  accepts  both  analog  and  digital  data  signals  from  various 
remote  sources  in  the  WU-2  aircraft  and  processes  them  into  a  digital  format  at  the 
rate  of  25  samples  per  second.  The  analog  data  signals  are  also  recorded  con¬ 
currently  in  analog  form  on  an  oscillogram.  Analog  signals  are  received  from  a 
gust  probe,  an  inertial  navigation  .ystem,  accelerometers,  thermometers,  altimeter, 
and  other  sensors.  The  gust  probe  provides  an  indication  of  airspeed  and  gust  forces. 
The  inertial  navigation  system  provides  information  on  aircraft  altitude,  acceleration, 
velocity,  and  horizontal  displacement. 

The  selection  of  clear  air  turbulence  samples  for  data  processing  starts  with  an  edit 
of  the  measurements  recorded  on  the  oscillogram.  Turbulence  samples  were 
selected  primarily  through  evaluation  of  the  records  of  acceleration  normal  to 
the  fuselage  axis.  This  acceleration  could  be  determined  to  within  0. 01  g.  Rapid 
variations  in  this  acceleration  for  ten  seconds  or  more,  with  amplitudes  exceeding 
+  0. 05  g  with  occasional  excursions  of  +  0. 10  g,  indicated  that  the  data  were  worth 
processing,  unless  records  indicated  that  the  aircraft  was  in  a  turn.  Small,  high- 
frequency  oscillations  in  the  gust  probe  records  normally  provided  further  evidence 
of  the  presence  of  turbulence.  Usually,  the  pilot  would  have  confirmed  the  presence 
of  turbulence  by  activating  the  oscillograph  event  marker. 

Each  edited  sample  was  characterized  by  a  subjective  description  in  words  in  addi¬ 
tion  to  the  estimated  eg  acceleration  peak  level.  In  general,  the  following  classifi¬ 
cation,  derived  originally  from  HICAT  pilot  comments,  was  observed: 
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+  0. 05  to  +  0. 10  g 
+  0.10  to +0.25  g 
+  0.25  to  +0.50  g 
+  0.50  to +0.75  g 
+  0. 75  or  greater 


Very  Light  (VL) 
Light  (L) 
Moderate  (M) 
Severe  (S) 
Extreme  (X) 


These  objective  criteria  were  not  complete  enough  to  eliminate  entirely  some  sub¬ 
jective  judgments  However,  the  limits  of  the  objective  criteria  were  sufficiently 
brood  that  the  role  of  subjective  judgment  in  the  selection  of  turbulence  cases  was 
relatively  minor.  For  Hie  purposes  of  this  report  the  turbulence  cases  presented 
by  Crooks,  Hoblit,  and  Prophet  (2)  were  accepted  without  modification. 

GRAPHS  AND  MAPS 

Once  the  specific  times  of  turbulence  were  selected,  records  from  all  the  various 
special  instruments  were  processed  by  techniques  given  by  Crooks  et  al.  (2 ). 

The  approximately  29  hours  of  turbulence  found  from  February  1965  to  March  1967 
required  processing  100  million  data  points.  The  end  products  of  this  data  pro¬ 
cessing  were  time  histories  and  power  spectral  densities.  In  addition,  flight  track 
maps  were  prepared. 

The  time  histories  were  presented  in  a  graphical  form  with  12  1/2  points  per  second 
along  the  time  axis.  Vertical,  lateral,  and  longitudinal  gust  velocities,  derived  gust 
velocity,  corrected  altitude,  ambient  air  temperature,  true  airspeed,  roll  angle, 
eg  longitudinal,  lateral  and  normal  acceleration  were  all  plotted  as  a  function  of  time. 
Including  selected  time  histories  from  the  newer  series  of  WU-2  flights  a  total  of 
322  time  histories  covering  approximately  1500  minutes  of  turbulence  were  available 
at  the  time  this  report  was  written.  A  sample  of  the  time  histories  is  given  in  Figure 
la  and  b.  .  Not  aU  the  time  histories  are  complete  because,  at  times,  one  or 
more  of  the  instruments  or  signal  recording  devices  did  not  operate  satisfactorily. 
Thus,  the  number  of  hours  of  prepared  time  histories  is  less  than  the  number  of 
hours  of  turbuience. 

Graphs  and  a  detailed  discussion  of  the  power  spectra  of  the  vertical,  lateral,  qnd 
longitudinal  gust  velocity  components  were  also  presented  by  Crooks,  Hoblit,  and 
Prophet  (2) .  ,  Some  of  the  graphs  of  the  power  spectra  for  the  latest  series  of 
WU-2  flights  were  also  available  for  use  in  the  present  study.  Figure  2  is  a  sample 
of  the  power  spectra  used  in  this  study.  The  problems  associated  with  classifying 
the  spectra  and  determining  relationships  between  spectral  types  and  meteorological 
and  terrain  conditions  are  discussed  later  in  this  secti  on. 

Flight  track  maps  for  the  present  study  were  prepared  by  two  methods.  In  some 
cases  they  were  prepared  directly  through  interpretation  of  the  records  obtained 
from  the  LN-3  navigational  system.  A  description  of  the  instruments  and  the  method 
of  interpreting  the  records  of  the  instrument  signals  is  given  by  Crooks,  Hoblit, 
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anfl  Prophet  (2).  Flight  track  reap*  were  based  upon  pile*  reports  when  the  LN~3 
equipment  did  not  operate  as  planned,  or  its  data  seemed  unreliable.  Samples  of 
flight  trade:  maps  are  presented  in  Figure  3i\  and  b. 

REVISED  TURBULENCE  MODEL 

Analysis  of  turbulence  measurements  requires  a  realistic  mathematical  model 
for  the  turbulence.  For  the  present  work*  the  model  used  for  the  previous 
study*  Hildreth,  Court*  and  Abrahms*  (4)  was  revised  alter  analysis  of  many 
such  models. 

Several  decades  of  work  on  theoretical  interpretation  of  experimental  wind  tunnel 
data  led,  by  the  late  1960’*,  to  the  general  adoption  in  aeronautical  engineering  of 
aa  expression  for  Hie  ^ectrum  of  vertical  gusts: 

*y(°>  =  (\  UH)  [l  +  3  (flL)2]  [l  +  (ftL)2]  “2  (1) 

o 

Here  $V(Q)  is  die  power  spectral  density  at  frequency  0  (radians  per  foot),  o 

Is  the  variance  of  the  vertical  motion,  or  mean  square  gust  velocity,  and  L  is  the 
longitudinal  scale  of  turbulence,  variously  interpreted  as  die  horizontal  distance  at 
which  die  correlation  between  horizontal  speeds  vanishes,  or  the  "effective  gust  size'.' 

This  turbulence  model  assumed  die  turbulence  .o  be  isotropic  and  the  correlation 
between  die  longitudinal  wind  components  (in  the  direction  of  flow)  to  decrease 
exponentially  with  separation.  It  was  used  in  many  studies  of  aircraft  response 
to  atmospheric  turbulence,  and  was  adopted  by  Hildreth,  Court,  and  Abrahms  (4) 
for  die  analysis  of  high- altitude  turbulence  spectra  obtained  by  Lockheed  in  1963-64, 
in  the  immediate  precursor  of  the  studies  reported  here. 

However,  Ibis  isotropic-exponential  model  has  certain  theoretical  drawbacks.  Most 
important  for  its  application  to  aircraft  design,  it  provides  that  the  spectral  power 
at  high  frequencies  varies  inversely  according  to  the  square  of  the  frequency.  But 
at  these  frequencies  the  well-established  similarity  theory  requires  that  the  variation 
be  according  to  die  -5/3,  not  -2  power.  Hence  von  Karin  an  ( 6 )  abandoned  the 
.  exponential  decrease  of  correlation  with  separation.  ’Tor  the  time  being”  he  pro¬ 
posed  ’to  interpolation  formula”  for  isotropic  turbulence  which  would  satisfy  the 
similarity  hypothecis. 

As  developed  in  detail  by  Houbolt,  Steiner,  and  Pratt  ( 7)  in  Appendix  B  of  their 
significant  report,  von  Karman’s  'Interpolation  formula”  leads  to 


(a2  L/n)  jl  +  8/3  (1.339  C1L)2]  [l  +  (1.339  CL)2] 


-11/6 


(2) 
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Figure  3a..  Flight  Track  Map  Over  Hawaii,  Flight  76. 


Figure  3b.  Flight  Track  Map  Over  Australia.  Flight  102 


Thus  the  final  term  has  an  exponent  of  -11/6  instead  of  -2  ,  the  coefficient  in 
the  second  term  is  8/3  instead  of  3 ,  and  a  numerical  constant,  involving  gamma 
functions,  is  introduced: 

r(f)  r(|)/r(|)  =  0.339... 

Extensive  comparisons  of  these  two  formulas  for  the  f>ower  spectra  density  of  vertical 
gusts,  by  Rhyne  and  Steiner  (8) ,  and  by  others,  led  Houbolk  Steiner,  and  Pratt 
to  conclude  that  the  von  Karman  representation  "best  fits  these  data". 

Consequently,  this  revised  model  has  also  been  adopted  for  the  analysis  of  atmo¬ 
spheric  turbulence  between  50, 000  and  80, 000  feet,  and  for  the  data  obtained  in  the 
experimental  flight  program.  However,  the  proper  value  of  the  scale  length,  L, 
in  Eq.  (2)  requires  investigation. 

Scale  Length 

Of  greater  interest,  both  theoretically  and  practically,  than  the  absolute  power  density 
spectrum  of  Eq.  (1)  and  (2)  is  the  relative  density,  §w(fl)  0^  »  obtained  by  dividing 

both  sides  by  the  total  variance  of  vertical  gusts.  This  relative  spectral  density 
depends  only  on  the  longitudinal  scale  length,  L  ,  defined  as 


<ux\+k>/(ul>  *  (3) 


the  integral  over  horizontal  distance  (x)  of  the  correlation  between  simultaneous 
horizontal  velocities  at  points  separated  horizontally  by  a  distance  k  .  This  correla¬ 
tion  decreases,  more  or  less  exponentially,  with  the  separation  k  ,  and  becomes 
essentially  zero  at  separations  greater  than  some  value  k_  .  Hence  the  integral  in 
Eq.  (3)  need  be  only  from  0  to  k^  ,  rather  than  to  ®  . 

*  * 

However,  the  correlation  can  be  computed  only  up  to  a  maximum  separation  k  , 
which  must  be  much  less  than  the  total  distance  over  which  measurements  are 
available,  and  usually  k*  «  k.  .  Hence  the  integral  cannot  be  computed  from  0 

u  * 

to  kft,  but  instead  must  be  truncated  at  k  ,  for  which  the  corresponding  lowest 
observable  frequency  is  Q  . 

In  practice,  L  is  evaluated  not  from  the  correlation,  but  by  solving  the  spectral 
density  function,  Eq.  (2),  to  yield,  Rhyne  and  Steiner,  ( 8 ) 

L  =  0.692  (av /a*)3  ^~3/3 -  ^  (4) 
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where  is  the  variance  of  vertical  motions  at  all  frequencies  and  o  *  the  vari¬ 
ance  over  the  observable  region,  assuming  that  almost  all  the  gust  energy  is  in  the 
observed  region,  end  0^  and  fi2  are»  respectively,  the  low  and  high  frequency 

ends  of  the  truncated  spectrum. 

In  addition  to  die  interpretations  as  maximum  distance  over  which  correlation  is 
appreciable,  and  as  maximum  effective  gust  size,  the  scale  length  L  has  another 
application  in  practice.  Logarithmic  differentiation  of  the  expression  for  the  rela¬ 
tive  spectral  density  (details  are  given  in  Appendix  1)  shows  that  the  spectral  density 
beecESi  constant  with  respect  to  0  at 


LQ  =  (1/2.187)0  =  0.0728X  (5) 


where  X  -  2tt/0  is  the  wavelength,  in  feet.  Thus,  the  scale  length  L  should  be 
about  7.3  percent  of  the  wavelength  at  which  the  spectrum  "levels  off".  (In  contrast, 
the  isotropic-exponential  model  used  previously,  Eq.  (1),  yields  a  criterion  of 
L_  *  1/1.7S20  *  0.0919X.)  Hence,  on  plots  of  relative  spectral  density  against 
frequency  (or  wavelength), values  of  Lmay  be  estimated,  although  some  extrapolation 
may  be  required.  Such  plots,  for  various  values  of  -L  ,  are  shown  in  Figure  4. 

Rhyne  and  Steiner  (8)  found  values  of  L  from  2, 000  to  "as  large  as  5, 000  feet 
for  severe- storm  turbulence  above  25, 000  feet  in  altitude'!  Reexamination  of  some 
clear-air  turbulence  data  by  Houbolt  et  al.  (9)  '^yielded  values  of  L  of  approxi¬ 
mately  3, 000  to  6, 000  feet  for  the  different  runs  and  methods  of  evaluation".  Hence 
they  concluded  "that  L  is  fairly  large,  apparently  of  the  order  of  5, 000  feet". 

However,  Hcblit  et  al.  (9)#  from  study  of  responses  of  various  subsonic  aircraft, 
accepted  2, 500  feet  as  the  best  estimate  of  scale  length  for  stratospheric  turbulence. 
Summarizing  this  work,  Stauffer,  Hildreth,  and  Crooks  (10)  pointed  out  that, 
with  realistic  values  for  root-mean- square  gusts  and  for  aircraft  operation,  the 
maximum  RMS  accelerations  for  supersonic  flight  increased  markedly  between 
assumed  scale  lengths  of  1, 000  and  2, 500  feet,  but  only  negligibly  when  the  length 
was  further  increased  to  5, 000  feet.  "SST  design,  based  on  a  turbulence  scale 
of  2, 500  feet, "  they  concluded,  "willnot  be  significantly  affected  if  additional 
research  indicates  a  higher  scale  to  be  in  order. " 

Reanalysis  of  data  from  the  first  HICAT  study,  Hildreth  and  Crooks  (11),  yielded 
spectral  curves  within  the  envelope  of  384  turbulence  spectra  from  B-66  flights  at 
low  levels,  and  overlapping  the  envelope  of  67  spectra  from  67  B-66  storm  penetra¬ 
tions,  as  well  as  the  envelope  of  10  TOPCA?  spectra,  Reiter  and  Burns  (12). 

In  aU  the  reanalyzed  HICAT  spectra,  the  scale  of  turbulence  apparently  was  slightly 
more  than  2, 500  feet. 


Reduced  Frequency  ft  (Radians/ft) 


INVERSE  WAVELENGTH,^-  (CYCLES/FT) 

Figure  4.  Relative  Spectral  Density  as  Function  ot  Longitudinal  Scale 
of  Turbulence  (L) 


Dynamic  Considerations 


The  scale  length,  L ,  also  is  related  to  other  atmospheric  properties.  Perturba¬ 
tion  theory  indicates  that  the  large  scale  static  stability,  wind  shear,  and  the  mean 
wind  speed  of  an  atmospheric  layer  are  the  important  parameters  that  determine 
the  character  of  wave  perturbations  (undulance)  in  that  layer  when  it  is  disturbed 
at  its  boundaries.  These  disturbances  may  result  from  low  level  air  flow  over  the 
terrain,  irregular  distribution  of  heat  sources,  (convection  in  particular),  as  well 
as.distuxbances  arising  from  dynamic  instability.  In  general,  the  short  waves  are 
external  waves  which  damp  out  quickly  with  height,  have  predominant  vertical  motion, 
and  may  be  unstable  if  the  wind  profile  is  right.  Long  waves  are  internal  waves  which 
are  generally  stable,  have  predominantly  horizontal  motion,  and  propagate  energy 
in  the  vertical. 

The  critical-wavelength  region  separating  the  short  external  waves  from  the  longer, 
internal- waves  appears  to  be  determined  principally  by  the  wind  speed  near  the 
boundary  and  the  static  stability  of  the  layer.  For  a  given  wind  speed,  the  greater 
the  static  stability  the  shorter  is  the  critical  wavelength.  For  a  given  static  stability 
the  greater  the  wind  speed  the  longer  the  critical  wavelength.  Disturbances  with 
frequencies  or  wavelengths  near  critical  values  may  cause  development  of  one  or 
more  resonant  waves,  with  large  vertical  amplitudes,  extending  through  deep  layers. 

Generally  the  static  stability  and  the  wind  combine  to  make  the  critical  wavelength 
in  the  stratosphere  on  the  order  of  1km.  However,  near  the  cold  side  of  a  jet  stream, 
where  the  static  stability  is  generally  less  but  the  wind  speed  still  high,  the  critical 
wavelength  may  be  as  great  as  a  100  km.  To  state  the  case  differently,  the  greater 
the  critical  wavelength  the  greater  the  probability  of  a  deep  layer  of  rough  air  with 
a  long  scale  length  and  greater  rms  gust  intensity. 

On  the  other  hand,  in  the  troposphere  the  critical  wavelength  generally  is  larger 
(10  -  50  km)  than  in  the  stratosphere.  Hence  any  resonant  mountain  waves  affecting 
the  stratosphere  would  most  likely  be  confined  to  thin  layers,  and  cause  horizontal 
internal  waves  to  develop,  except  possibly  over  and  on  the  cold  side  of  jet  streams, 
as  indicated  above.  The  longer  external  waves  would  also  be  confined  to  thin  strato¬ 
spheric  layers  but  could  cause  internal  waves  to  develop.  The  shorter  external 
waves  in  the  troposphere  could  excite  external  waves  in  the  stratosphere,  but  both 
data  and  theory  indicate  that  the  energy  in  these  tropospheric  external  waves  would 
probably  be  negligible  at  the  tropopause  except  when  the  mountains  are  high,  the 
tropopause  low  and  the  wind  speed  large. 

Disturbances  with  wavelengths  longer  than  those  critical  in  the  troposphere  should 
propagate  energy  upward.  Theoretical  estimates,  Eliassen  and  Palm  (13),  indicate 
that  waves  shorter  than  about  20  km  (usually  external  waves)  are  effectively  reflected 
by  the  upper  troposphere.  Short  internal  waves  ( <  50  km)  have  approximately  up 
to  three  fourths  of  the  energy  reflected  while  longer  waves  have  progressively  less 
energy  reflected  in  the  upper  troposphere. 

Disturbances  from  low  level  heat  sources  would  have  similar  scale  effects  on  the 
stratosphere.  The  importance  of  convection  in  perturbing  stable  layers  has  been 
investigated  by  Townsend  (14).  As  the  scale  of  thermals  affecting  the  stratosphere 


is  on  the  order  of  1  to  10  km,  external  waves  and  resonant  waves  could  be  excited 
as  well  as  internal  waves. 

These  considerations  show  that  the  important  meteorological  parametTf!  are: 

(a)  The  critical  wavelength  in  the  stratosphere 

(b)  The  tropospheric  critical  wavelength  for  terrain  effects,  and 

(c)  A  measure  of  convective  instability  (such  as  the  Showalter  Index). 

The  geophysical  parameters  include  a  measure  of  the  terrain  roughness,  as  well  as 
a  height-to-width  ratio  for  determining  resonant  waves.  For  heating  effects,  a 
measure  of  the  irregularity  of  the  surface  vegetation  character  should  be  important. 
All  of  these  eventually  may  provide  a  direct  means  for  estimating  the  scale  length, 

L  ,  and  other  characteristics  of  a  high  altitude  turbulence  model.  But  until  they 
have  been  investigated  in  detail,  the  most  useful  approach  to  analysis  of  high  altitude 
spectra  appears  to  be  the  application  of  the  von  Karman  model,  Eq.  (2),  assuming 
L  *  2, 500  feet. 

CLASSIFICATION  OF  SPECTRA  OF  GUST  VELOCITIES 

In  the  inertial  subrange  the  power  spectral  density  should  depend  on  the  frequency 
Q  and  the  energy  dissipation  rate  e  according  to 


=  a  e2/3ffp  (6) 


where  a  is  a  universal  constant  and  p  is  the  slope  of  the  power  spectral  density 
curve.  Several  authors,  e.g.,  Reiter  and  Burns  (15);Pinus,  Reiter,  Shur,  and 
Vinnichenko  (16),  have  obtained  an  experimental  value  of  p  =  -5/3  in  the  inertial 
subrange  in  agreement  with  that  derived  theoretically  by  Komogorov  ( 17) .  Other 
values  for  p  are  of  course  possible,  such  as  -2,  Eq.  (1),  depending  on  the  assump¬ 
tions  made  in  deriving  an  expression  for  $(fl)  .  The  HICAT  spectra  show  significant 
deviations  from  both  the  -5/3  and  -2  laws. 

Spectral  Slopes 

Energy  spectra  derived  from  time  histories  of  gust  velocities  are  commonly  pre¬ 
sented  by  plotting  the  power  spectral  density  divided  by  the  inverse  wavelength  in 
units  of  (ft/seclVcycles/ft  as  a  function  of  the  inverse  wavelength  in  cycles/ft 
(Figure  2).  However,  the  spectra  available  from  HICAT  flights  54  through  265  have 
been  normalized  with  respect  to  the  energy  by  dividing  the  ordinate  by  the  variance 
of  the  vertical  gust  velocity  for  wavelengths  less  than  2000  ft.  The  ordinate  is  hence 
no  longer  the  same  measure  of  turbulent  energy  as  before  and  is  in  units  of  (cycles/ftf. 
The  normalization  process  has  effectively  reduced  to  a  minimum  the  differences  in 
energy  for  X  <  2000  ft  and  facilitated  the  comparison  of  spectral  shapes,  especially 
for  X  >  2000  ft.  The  2000  ft  cutoff  was  selected  in  order  that  all  available  spectra 
could  be  included. 
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Following  normalization,  spectra  for  the  vertical  component  of  the  gust  velocity 
(Oy)  were  classified  into  several  groupings  according  to  shape.  Only  those  cases 
were  considered  where  RMS  >  0.63  because  pilot-induced  control  move¬ 
ments  or  instrument  malfunctions  affect  the  spectra  shapes  appreciably  at  lower 
variances. 

Figures  5  through  7  present  the  Uy  spectra  for  various  categories  number  A  through 
F*.  The  p  *  -5/3  and  p  *  -2  lines  have  been  entered  for  comparison  purposes. 
The  solid  curved  lines  are  composites  of  the  several  cases. 


Ths  distribution  of  spectra  by  classes  is: 

CLASS  A  B  C  D 

NUMBER  21  8  13  15 


Unclassified 

9 


In  general,  the  average  slope  of  the  spectra  increases  with  class  from  A  to  F*. 

There  are,  however,  singularities  in  the  spectral  shapes  which  will  be  described 
in  more  detail. 

The  spectrum  curve  for  classes  A  and  D  has  a  value  for  p  much  less  than.  -5/3  for 
X  <  1000  and  2000  ft,  respectively.  Class  C  spectra,  in  addition,  has  p  <  -5/3  in 
the  approximate  range  1000  ft  >  X  >  150  ft;  E  and  F  spectral  slopes  are  less  than 
-5/3  for  X  <  2000  ft.  Approximately  80%  of  the  classified  spectra  deviate  towards 
shallowed  slopes  than  -5/3  for  X  <  1000  ft,  contrary’to  the  results  of  others  (Pinus, 
Reiter,  Shur,  and  Vinnichenko,  (16) ;  Reiter  and  Burns,  (13),  for  flights  at  commer¬ 
cial  jet  altitudes  (i.e. ,  below  40, 000  ft). 

Approximately  70%  of  the  spectra,  including  all  the  classes  with  the  exception  of  D, 

F,  and  F*,  show  a  tendency  towards  constant  or  decreasing  p  with  increasing  X . 

B  spectra  especially  show  a  flattening  out  at  longer  wavelengths.  Some  portions  of 
the  E*,  F,  and  F*.  curves,  43%  of  the  total  classified  spectra,  possess  p  >  2  for 
X  >  900  ft.  The  long  wave  portions  of  the  spectra  for  runs  102  (16),  114  (14),  and 
219  (3)  in  the  F*  category  have  -3  slopes  which  agrees  with  the  findings  of  Shur  (1962). 

A  total  of  9  individual  cases  vary  substantially  from  the  categories  A  through  F*  and 
are  grouped  separately  (Figure  8).  Spectra  from  flight  100,  run  16  and  flight  257,  run 
run  13  (Figure  8a)  fall  within  the  C  class  for  X  >  800  ft  but  deviate  from  this  at 
small  X  .  Flight  258,  run  8,  10  and  flight  264,  run  14  (Figure  8b)  are  characteristic 
of  the  B  class  at  shorter  X  and  E  or  E*  at  X  >  900  ft.  The  remaining  four  cases 
(Figure  8c)  do  not  appear  to  fall  into  any  basic  category  or  combination  of  categories. 

The  spectra  illustrated  in  Figure  9a  provide  an  example  that  indicates  that  flight 
orientation  with  respect  to  the  wind  has  no  significant  effect  upon  the  shape  of  spectral 
curve.  The  runs  were  only  6  minutes  apart  and  over  the  same  course  exfeept  run 
107  j(9)  was  into  the  wind  and  107  (10)  with  the  wind.  The  wind  speed  was  around 
65  KTS  at  the  time. 

Figure  9b  shows  how  there  can  be  but  small  variations  in  the  spectral  shapes  for 
several  oases  over  a  wide  range  of  distance.  Eleven  samples  were  taken  in  a  600 
mile  path  along  the  lee  side  of  the  Great  Dividing  Range  in  Southeastern  Australia. 
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The  lateral  (U^)  and  longitudinal  (Up)  gust  velocity  components  show  a  tendency 
towards  less  variation  in  their  spectral  shapes  and  are  not  considered  in  detail. 
They  appear  to  possess,  on  the  average,  slightly  steeper  slopes  than  the  Uy 
spectra.  A  comparison  between  the  three  components  is  shown  for  two  flights  in 
Figures  lOa-b.  Flight  76  is  over  water,  near  the  Hawaiian  Islands,  flight  114  is 
on  the  lee  side  cf  the  Sierra  Nevada  range  in  Califomia-Nevada.  The  spectra  agree 
in  shape  reasonably  well  for  both  cases,  following  the  -5/3  law. 

Terrain 

From  a  total  of  126  Uy  spectra,  15%  were  over  water,  30%  over  relatively  flat 
plains,  and  55%  over  or  reasonably  close  to  mountainous  terrain.  The  distribution 
of  classes  by  terrain  (Table  It)  shows 

TABLE  H 

NUMBER  OF  CASES  FOR  EACH  Uy  SPECTRA  AND  TERRAIN 
CLASSIFICATION  FDR  FLIGHTS  54-265  (126  CASES) 


Class 

Water 

Flatland 

Mountain 


Unclassified 

1 

3 

5 


Total 

19 

38 

69 


Grouping  classes  A  through  D  into  one  category  and  E  through  F*  into  a  second  cate¬ 
gory,  the  distribution  by  terrain  shows  86%  of  the  flatland  spectra  in  the  A  through 
D  class  and  72%  of  the  mountain  spectra  in  the  E  through  F*  class.  The  water  spectra 
are  divided  equally  between  the  two  classes.  However,  the  sample  size  in  this  case 
is  small. 

Meteorological  Conditions  Associated  with  Spectral  Classes 


The  analysis  of  the  meteorological  conditions  associated  with  the  various  spectral 
classes  is  limited  to  flights  54  to  164  because  of  incomplete  data  analysis  for  later 
flights.  For  these  flights ,  investigation  indicated  that  there  was  no  significant  rela¬ 
tionship  between  the  spectral  classes  and  various  categories  of  stability,  wind  shear 
jet  strength,  Richardson  number,  and  proximity  to  thunderstorms.  No  conclusive 
evidence  could  be  drawn  from  the  data.  The  reasons  for  this  are  as  follows: 

1.  The  original  sample  size  is  too  small. 

2.  Radiosonde  and  rawinsonde  data  are  lacking  in  many  cases  at  the  heights 
flown. 

3.  It  is  difficult  to  resolve  scale  size  differences  between  the  meteorological 
variables  and  turbulence  samples. 

4.  There  appears  to  be  an  overlap  between  the  effect  of  terrain  and  of  various 
meteorological  conditions  on  the  spectral  shapes. 
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5.  Turbulence  samples  are  limited  in  their  geographical  extent.  Also,  the 
method  of  pattern  flying  has  resulted  in  several  turbulence  samples  for 
very  few  types  of  synoptic  situations. 

6.  An  individual  spectrum  represents  an  average  of  the  eddy  conditions  over 
the  length  of  the  run  and  its  shape  varies  depending  on  the  time  period 
selected  for  the  rim.  Consequently,  spectra  derived  frc~a  data  obtained 
within  an  area  associated  with  a  specific  synoptic  situation  can  exhibit 
diverse  shapes. 


SECTION  IV 


UPPER  AIR  DATA 


Two  principal  sources  of  meteorological  data  were  used  in  the  present  study*  One 
source  was  the  records  of  the  ambient  temperature  presented  in  the  time  histories 
for  the  WU-2  flights.  In  these  time  histories  the  ambient  temperature  is  shown  as 
a  function  of  time  by  12. 5  points  for  each  second  along  the  time  axis.  The  response 
ox  the  thermometer  was  such  that  temperatures  variations  slower  than  one  cycle 
per  second  could  be  followed.  The  second  source  of  meteorological  information 
was  the  twice  daily  (once  dally  in  Australia  and  New  Zealand)  temperature,  pressure 
and  wind  measurements  reported  by  the  various  national  meteorological  organizations. 
The  number  of  observations  made  per  unit  area  and  the  type  of  instruments  used 
varied  from  one  HICAT  search  area  to  another.  Specifically,  the  areas  around  the 
Hawaiian  Islands  and  around  the  Panama  Canal  Zone  have  -relatively  few  upper 
air  observations  at  altitudes  of  40*  000  to  VO,  000  ft.  In  general,  the  national 
meteorological  organizations  make  their  observations  of  the  upper  atmosphere  near 
G0Z  and  12Z  times.  The  time  and  space  interval,  between  the  routine  meteorological 
obae-rs  aliens  and  the  turbulence  encountered  by  the  WU-2  aircraft  varied  significantly. 
Only  rarely  was  high  altitude  turbulence  reported  over  a  meteorological  station  within 
an  hour  of  the  measurement  made  by  an  instrument  carried  aloft  by  a  balloon. 

Figures  11-17  illustrates  the  relative  numbers  and  locations  of  the  meteorological 
stations  at  which  temperatures,  pressures,  and  winds  are  measured  at  WU-2  fligit 
levels  for  each  of  the  HICAT  search  areas. 
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Figure  12.  Upper  Air  Station  Locations.  New  Zealand 
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SECTION  V 


DISCUSSION 


AN  INTERPRETATION  OF  THE  RADIOSONDE  TEMPERATURE  PROFILES  IN 
TERMS  OF  VERTICAL  WAVE  MOTION  IN  THE  ATMOSPHERE 

Radiosondes  are  carried  aloft  by  balloons  and  the  balloons  have  a  horizontal  com¬ 
ponent  of  motion  that  is  associated  with  the  wind.  Hence*  the  temperature-height 
profiles  obtained  from  radiosonde  data  are  actually  temperature  profiles  along  a 
three  dimensional  curve  that  isf  in  general*  not  a  vertical  line  extending  through 
the  point  of  origin.  For  many  purposes  it  suffices  to  assume  that  the  radiosonde 
record  provides  an  adequate  approximation  to  a  vertical  sounding.  In  the  following 
discussion  it  will  be  shown  that  under  certain  conditions,  the  departure  of  the  radio¬ 
sonde  path  from  the  vertical  may  result  in  a  temperature  profile  that  shows  signifi¬ 
cant  changes  in  lapse  rate  which  would  not  be  indicated  if  the  radiosonde  rose  ver¬ 
tically,  It  will  be  shown  that  these  changes  in  lapse  rate  may  be  interpreted  in 
terms  of  wave  motion  in  the  atmosphere.  Such  an  interpretation  of  the  radiosonde 
temperature  profiles  provides  a  physical  explanation  of  the  apparently  conflicting 
statements  found  in  the  Hterature  relative  to  the  correlation  between  high  altitude 
clear  air  turbulence  and  temperature  lapse  rate. 

The  analysis  that  follows  is  an  amplification  of  a  discussion  given  by  Scorer  (18). 
Briefly,  Scorer  considered  the  case  of  a  balloon-borne  radiosonde  rising  at  a  con¬ 
stant  rate  relative  to  the  air  and  also  moving  horizontally  through  a  stationary  wave 
in  the  atmosphere  (mountain  wave).  He  pointed  out  that  the  temperature  indicated 
by  the  radiosonde  sensor  at  any  given  height  would  be  a  function  of  position,  i.  e. , 
the  phase  angle  of  the  atmospheric  wave  at  that  particular  height.  The  radiosonde 
traveling  through  the  wave  would  indicate  the  existence  of  temperature  inversions 
that  would  not  be  present  for  a  sounding  that  had  followed  a  truly  vertical  path. 

This  may  be  illustrated  by  a  hypothetical  case  presented  in  Figures  18a  and  18b. 

In  this  case,  it  is  hypothesized  that  a  stationary  wave  motion  of  constant  amplitude 
and  wavelength  is  initiated  in  an  isothermal  layer.  Due  to  the  wave  motion,  the 
position  of  the  potential  temperature  surfaces  might  be  as  illustrated  in  Figure  18a. 

If  a  radiosonde  enters  the  wave  at  point  ,TA?5  and  follows  the  path  ABCD,  the  tempera¬ 
ture  profile  as  indicated  by  the  radiosonde  would  then  be  a  nearly  sinsoidal  curve 
which  in  a  teletype  transmission  would  be  approximated  by  the  straight  lines  given 
in  Figure  18b.  A  true  vertical  sounding  would  indicate  an  isothermal  atmosphere 
with  temperatures  varying  from  a  minimum  in  the  ridges  (maximum  adiabatic 
cooling)  to  a  maximum  in  the  troughs  (maximum  adiabatic  heating).  Of  course, 
if  the  wave  were  moving,  then  even  a  true  vertical  sounding  would  indicate  layers 
of  varying  lapse  rate. 

It  is  proposed  here  that  the  simple  stationary  wave  model  provides  a.physical  explan¬ 
ation  for  the  apparent  relationship  between  the  occurrence  of  high  altitude  clear  air 
turbulence  and  the  presence  of  temperature  inversions  in  the  radiosonde  temperature 
profiles  as  described  in  the  literature.  Further,  this  simple  model  may  be  used  to 
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Figure  18.  Hypothetical  Case  of  Balloon  Ascent  Through  a 
Wave  and  Resulting  Temperature  Sounding 
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explain  the  disagreements  on  the  relative  position  of  the  observed  turbulence  with 
respect  to  the  inversion.  For  example,  Haymond  (19)  claims  that  the  turbulence 
occurs  in  an  inversion  layer  and  others  claim  that  turbulence  occurs  between  inver¬ 
sion  layers.  The  analysis  outlined  above  indicates  that  the  height  of  an  inversion 
layer  on  the  radiosonde  temperature  profile  may  simply  be  a  function  of  the  balloon's 
position,  and  that  the  observed  clear  air  turbulence  may  be  the  result  of  wave  motion 
in  the  atmosphere. 

CALCULATION  OF  MEAN  SQUARE  VERTICAL  ACCELERATION  FROM  RADIO¬ 
SONDE  DATA 

The  mean  square  vertical  accelerations  associated  with  clear  air  turbulence  were 
calculated  starting  with  the  postulate: 

a.  The  wave  motion  in  a  stably  stratified  atmosphere  layer  has  constant 
phase  and  amplitude  with  altitude. 

b.  Only  dry  adiabatic  thermodynamic  processes  are  significant. 

c.  The  balloon  rises  with  a  constant  speed  of  5. 08  m/sec  (1000  ft/min) 
relative  to  the  air. 

d.  The  wind  is  normal  to  the  crest  2nd  trough  lines  and  for  a  stationary 
wave  is  equal  to  u,  the  mean  speed  in  the  layer. 

Under  these  conditions  the  horizontal  distance  that  the  radiosonde  will  travel  in  going 
from  a  minimum  to  a  maximum  temperature  (maximum  to  a  minimum)  is  given  by 


x  =  X/2  =  u  Az/wq 


(?) 


where  X 


u 

Az 


w 


o 


wavelength  of  the  wave  in  the  atmosphere 

horizontal  component  of  the  wind  relative  to  the  wave 

altitude  interval  between  observed  minimum  and  maximum  temperature 

rate  of  rise  of  balloon  with  respect  to  air 


For  a  sins  wave  the  vertical  displacement,  z  ,  may  be  expressed  by 


z  =  Ah  sill  2tt  (t/ t^) 


(8) 


where  Ah  =  amplitude  of  the  wave 
Tj  *  wave  period 

t  «  time 
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The  expression  for  the  vertical  velocity  is  obtained  by  differentiating  Eq.  (8)  with 
respect  to  time: 


dz/dt  -  w  =  2rr  ^  cos  2tt  ~ 
T1  T1 


(9) 


If  an  airplane  were  to  fly  horizontally  through  such  a  wave  it  would  be  subjected  to 
vertical  motions  of  the  atmosphere  given  by 


w 


2tTAh 

T1 


cos  2tt 


X 

t2 


(10) 


where  t2  Is  the  period  of  the  wave  as  experienced  by  the  aircraft  or 


t2 


(10a) 


and  v  is  the  speed  of  the  aircraft  relative  to  the  speed  of  the  wave. 

8 

The  acceleration  of  the  air  through  which  the  aircraft  is  flying  is  obtained  by 
differentiating  Eq.  (10) 


dw/dt  =  ••  sin  (2tt  ~) 

T1  t2  T2 

Squaring  both  sides  we  have 

(W 

T  T  2 

1  t2 


(11) 


(12) 


To  obtain  the  mean  square  acceleration  both  sides  of  Eq.  (12)  are  integrated  over 
the  limits  t  ®  0  and  t  *  T2/r  and  divided  by  Tg/4  yielding 
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(13) 


(dw/dt)2  =  ZZ-iQ.)2 

2  2 

Ti  '2 


and  taking  the  square  root  lo  obtain  the  root  mean  square  acceleration 


RMS  =  J  (dw/dt)2  = 

T1  t2 


(14) 


Substituting  *  X/u  and  t2  3  X/vg 


RMS  =  2  J2  tt  fr.,LvB 
X2 


(15) 


or  since  X  *  from  Eq.  (7) 


Ah  v„  w' 


RMS  vertical  acceleration  =  rr  /JIT  - - 

V  a  (Az)2 


(16) 


This  equation  implies  that  the  RMS  vertical  acceleration  of  the  atmosphere  with 
respect  to  the  aircraft  is  directly  proportional  to  the  speed  of  the  aircraft  relative 
to  the  speed  of  the  wave  and  to  the  amplitude  (Ah)  of  the  wave  but  inversely  pro¬ 
portional  to  the  square  of  the  depth  of  the  layer  between  the  observed  maximum 
and  minimum  temperature.  Hence,  with  other  things  being  equal,  a  sounding 
depicting  fluctuations  in  the  temperature  gradient  over  relatively  shallow  layers 
implies  more  severe  turbulence  than  when  the  same  temperature  changes  take 
place  through  deeper  layers. 

A  sample  computation  is  given  for  flight  test  114  (29  Aug.  1966).  On  this  day  exten¬ 
sive  light  to  moderate  turbulence  was  found  on  the  lee  side  of  the  Sierra  Nevada 
mountains  from  just  south  of  Reno  to  Las  Vegas.  The  Las  Vegas  radiosonde  obser¬ 
vation  for  30  Aug.  1966  0000Z  is  shown  for  in  Figure  19.  The  turbulence  was 
reported  between  53  to  57, 000  It,  corresponding  well  with  the  fluctuations  on  the 
radiosonde  temperature  profile.  A  stationary  wave  is  assumed.  The  following 
calculation  was  made  from  the  sounding  for  the  three  layers  A,  B,  and  C. 
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Ah  *  483  ft  (147  m) 

Ax  «  1840  ft  (564  m) 
u  *  25  knots  (12.5  m/sec) 
vg  »  400  knots  (206  oi/sec) 

X  «  s  9106  ft  ^776  m) 


RMS  vertical  acceleration 


2  /T  ^(147)02.5^06) 
(2776)2 


1.37  m/sec2 


Assuming  the  aircraft  responds  to  all  vertical  Sir  motion,  then  this  RMS  vertical 
acceleration  can  be  compared  with  the  observed  RMS  incremental  eg  normal  acceler¬ 
ation,  Aaj- ,  of  the  aircraft  which,  when  expressed  in  terms  of  the  acceleration  of 
gravity,  is  computed  to  be 


RMS 


1.37 

9.8 


0.14 


The  observed  RMS  Aa^/g  mean  values,  RMS  AG  in  Crooks,  Hoblit  and  Prophet  (2) 
near  Las  Vegas  ranged  from  .  051  to  .  078. 

The  aircraft  did  not  fly  near  any  other  radiosonde  stations  on  this  day,  but  on  the 
following  day.  Test  115  30  Aug.  1966,  in  addition  to  finding  turbulence  over  the 
Sierra  Nevada  Mountains,  tight  turbulence  was  also  found  at  61  thousand  feet  over 
Ely.  The  Ely  sounding  nearest  this  time  is  shown  in  Figure  20.  However,  no  tur¬ 
bulence  was  found  from  just  west  of  Ely  to  Reno.  The  nearest  radiosonde  to  this 
no  turbulence  area  is  Winnemucca  shown  in  Figure  21.  Note  the  contrast  compared 
to  the  Ely  sounding. 

Sample  calculations  for  other  areas  have  also  been  made.  On  16  May  1966  (Test  76) 
widespread  light  and  moderate  turbulence  occurred  around  the  large  island  of 
Hawaii  at  about  58, 000  ft.  'The  Hilo  sounding  at  this  time  is  shown  in  Figure  22 
and  according  to  the  proposed  theory,  indicates  that  turbulence  is  likely  from 
50, 000  ft  to  just  above  60, 000  ft.  The  RMS  AaN/g  value  calculated  from  the 
sounding  is  .  075.  The  observed  mean  values  for  all  runs  in  the  Hilo  area  ranged 

f-rvvm  fA  AK1  c* 


On  15  Apr  1966  (Test  58),  no  turbulence  could  be  found  anywhere  in  the  Hawaiian 
area  over  the  altitude  range  50  to  68, 000  ft.  The  Hilo  and  Lihue  soundings  at  this 
time  are  shown  in  Figures  23  and  24,  respectively.  As  may  be  noted,  no  wave 
motion  is  indicated  at  any  altitude  from  50  to  68, 000  ft. 


On  28  Jul  1966  (Test  102),  moderate  turbulence  was  quite  widespread  to  the  lee  of 
the  great  dividing  range  in  Australia  while  no  turbulence  was  encountered  t  >,  the  west 
of  the  range.  An  extensive  area  of  moderate  turbulence  was  found  at  60, 000  ft  over 


I 


45 


Brisbane.  The  radiosonde  at  this  time,  shown  in  Figure  25  indicates  probable  wave 
motion  over  the  interval  56  to  65,  000  ft.  for  the  layer  58.5  to  62.5  thousand  feet, 
the  computed  RMS  Aa^/g  value  is  .  115.  '.'he  observed  mean  values  for  all  runs 
in  the  Brisbane  area  range  from  .  038  to  .  095, 

On  the  same  day  an  area  of  moderate  turbulence  was  found  at  about  63,  000  ft  just 
north  of  Williamstown  and  patches  of  lighter  turbulence  at  65  to  66,  000  ft  to  the 
south.  The  Williamstown  sounding,  shown  in  Figure  26  indicates  probable  wave 
motion  irom  57  to  66,  000  ft.  The  RMS  Aa^/g  value  calculated  from  the  sounding 
for  the  layer  62  to  66,  000  ft  was  .  082.  The  observed  mean  values  for  all  runs  in 
the  Williamstown  area  ranged  from  .  050  to  .  059. 

No  turbulence  was  found  to  the  west  of  the  range.  The  soundings  for  Cobar  and 
Wagga,  shown  in  Figui'es  27  and  28  respectively,  indicate  no  wave  motion  at  flight 
altitudes,  altnough  moderate  turbulence  was  found  over  the  mountains  less  than 
100  miles  east  of  Wagga. 

It  appears  from  these  exan  pies,  summarized  in  Table  III,  pg.  64,  that  when  the  radio¬ 
sonde  observations  are  taken  fairly  close  to  tM  time  of  flight  good  estimates  ci  mean 
turbulence  intensity  can  be  obtained  from  the  meteorological  data.  In  the  few  rare 
instaices,  notably  in  Australia,  when  special  radiosonde  observations  were  taken 
six  hours  ap art,  it  is  evident  that  the  temperature  structure  can  change  quite 
drastically.  Hence,  in  areas  where  only  the  regular  12-hourly  observations  are 
made  and  the  flight  time  is  midway  between,  lower  correlations  between  observed 
and  calculated  values  can  be  expected.  The  distance  to  which  a  particular  radio¬ 
sonde  might  be  applicable  will  undoubtedly  be  greater  in  areas  of  uniform  terrain, 
whether  rough  or  smooth,  as  compared  to  an  area  in  which  there  is  an  abrupt 
change  in  terrain  profile,  as  was  noted  in  Eastern  Australia. 

The  assumptions  used  in  deriving  this  simple  wave  model  should  act  to  maximize 
the  computed  values  for  RMS  Aa^/g  .  For  instance,  if  the  wave  train  is  not 
stationary  but  moving  downwind  at  some  speed  less  than  the  mean  wind,  the  vertical 
velocities  and  accelerations  would  be  smaller.  If  the  depth  of  the  layer  or  the 
relative  wind  speeds  are  so  small  that  the  balloon  does  not  have  time  to  travel  down¬ 
wind  relative  to  the  wave  a  distance  equal  to  half  a  wavelength  while  passing  up 
through  the  wave  layer,  the  computed  values  for  Ah  and  Az  in  Eq,  (16)  will  be 
too  small,  and  since  Az  is  squared  and  in  the  denominator,  the  computed  RMS 
Aan/g  would  again  be  on  the  large  side.  Moreover,  it  is  assumed  that  the  aircraft 
responds  to  all  vertical  air  motion,  i.e.,  there  is  no  vertical  motion  of  the  air 
relative  to  the  aircraft.  Actually,  due  to  inertia,  the  aircraft  response  is  usually 
less. 

A  MODEL  FOR  CIEAR  AIR  TURBULENCE  BASED  ON  SLOPING  ISENTROPIC 
SURFACES 

The  model  for  identifying  turbulent  layers  and  computing  approximate  values  for 
root  mean  square  vertical  accelerations  presented  in  the  previous  section,  is  appli¬ 
cable  only  to  situations  when  the  turbulence  happens  to  be  occurring  at  the  time  and 
place  of  radiosonde  ascent.  Turbulence  will  sometimes  occur  within  *he  relatively 
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wic<  ly  spaced  radiosonde  network  but  not  necessarily  near  enough  to  any  station 
to  be  analyzed  from  the  sounding.  Therefore,  there  appeared  to  be  a  need  to 
develop  a  model  for  identifying  potential  turbulence  from  the  analyzed  fields  of 
temperature  and  pressure  over  relatively  large  areas  as  depicted  on  upper  air 
synoptic  charts. 

It  is  assumed  that  turoulence  in  the  stratosphere  is  associated  with  vertically  rising 
and  descending  air  and,  due  to  the  stable  stratification,  any  vertical  motion  must 
be  forced,  i.e.,  energy  consuming.  Examples  of  mechanisms  which  could  cause 
forced  vertical  motion  are  convection,  mountains,  fronts,  and  mass  divergence 
associated  with  synoptic  pressure  systems.  In  the  case  of  convection,  we  are 
concerned  with  the  so-called  negative  buoyancy  regior  above  active  convective 
cells.  The  role  of  frontal  systems  and  mountains  are  similar  in  that  both  obstruct 
the  air  flow  and  may  cause  it  to  rise  on  the  windward  and  descend  on  the  leeward 
side.  Synoptic  pressure  systems  may  involve  large-scale  horizontal  mass  diver¬ 
gence  that  is  associated  with  vertical  adjustment  of  air  masses,  and  the  more  intense 
systems  extend  upward  into  the  stratosphere.  Continuity  considerations  require 
that  diverging  or  converging  horizontal  wind  fields  be  associated  with  vertical 
motion,  ns  the  jet  stream  itself  may  be  subjected  to  vertical  displacements  as 
it  passes  around  a  trough  or  ridge  aloft,  the  jet  appear  to  play  no  direct  role  in 
producing  vertical  motion. 

The  proposed  model  is  one  in  which  an  initial  perturbation  imparts  a  vertical  com¬ 
ponent  to  the  air  motion.  This  perturbation  is  presupied  to  be  produced  by  some 
intermediate  to  large-scale  phenomena  ranging  from  a  few  miles  in  diameter  in 
the  case  of  forced  convective  activity  to  100  miles  or  more  in  the  case  of  synoptic 
scale  phenomena.  As  adiabatic  temperature  changes  are  assumed,  then  perturba¬ 
tions  of  vertical  displacement  should  be  associated  with  perturbations  in  the  tempera¬ 
ture  field  along  the  direction  of  flow.  Methods  of  analyzing  these  perturbations  can 
therefore  be  based  on  the  pattern  of  isotherms  on  a  constant  pressure  surface  or 
from  a  height  contour  analysis  of  an  isentropic  surface. 

As  a  consequence  of  the  vertical  displacement  of  air  masses  affected  by  the  initial 
perturbations,  buoyance  and  Inertial  forces  produce  secondary  perturbations  or 
waves.  The  period  of  these  secondary  perturbations  is  a  function  of  the  vertical 
stability  and  also  vertical  wind  shear.  The  resulting  wavelengths,  being  cf  the 
order  of  a  few  miles,  are  short  enough  to  produce  significant  vertical  accelerations 
on  aircraft  flying  through.  In  the  ~vent  the  period  of  the  initial  perturbation  corres¬ 
ponds  to  wavelengths  of  this  size,  as  is  likely  in  convective  activity,  then  the  initial 
perturbation  itself  can  cause  vertical  accelerations.  However,  for  the  larger  scale 
phenomena,  the  wavelengths  would  he  much  too  long  to  cause  significant  vertical 
acceleration  on  aircraft.  It  should  be  noted  that  according  to  this  model  the  energy 
does  not  eascade  down  from  the  large  scale  initial  perturbations  to  the  smaller 
scale  secondary  perturbations  and  thereby  yield  a  continuous  spectrum  of  wave¬ 
lengths,  On  the  contrary,  it  is  hypothesized  that  the  secondary  perturbations  wldch 
contain  the  maximum  vertical  component  of  turbulent  energy  occur  at  discrete  fre¬ 
quencies  depending  upon  trie  vertical  stability  and  wind  shear. 
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The  Model 


In  a  stably  stratified  layer,  air  which  is  forcibly  displaced  in  the  vertical  will  be 
subjected  to  deceleration  due  to  buoyance  forces.  The  air  eventually  stops  rising 
or  decending  and  is  accelerated  back  toward  the  level  of  origin  or  equilibrium  level. 
However,  inertia  carries  the  air  past  this  level  and,  assuming  no  loss  of  energy, 
the  air  oscillates  about  the  equilibrium  level  with  amplitude  h  and  frequency  v  . 

The  frequency  is  related  to  the  period  t  ,  defined  as  the  time  required  for  one 
complete  oscillation,  by 


V  -  2  Ti  /' 


(17) 


and  the  period  is  related  to  the  mean  wind  speed  and  wavelength  by, 


T  =  X/U 


(18) 


Assuming  only  adiabatic  temperature  changes,  the  amplitude  h  will  equal  the  max¬ 
imum  vertical  displacement  of  a  potential  temperature  surface  from  its  equilibrium 
level.  The  vertical  displacement  z  in  a  sine  wave  is  related  to  amplitude  and  time 
by 


z 


h  sin 


2nt 

T 


(19) 


A r  expression  for  vertical  velocity  is  obtained  by  differentiating  Eq.  (19)  with  respect 
tc  time. 


dz/dt  =  w  =  ~  h  cos 

T  T 


(20) 


from  which 


w  = 


4tt 


,2  2  2  TTt 

h  cos  - - 


(21) 
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Integrating  both  sides  of  Eq.  (21)  with  respect  to  time 


t  =  -/4 

[  2  4 

J  W  dt  =  - 
t  =  0 


and  dividing  both  sides  by  ~/A  ,  the  expression  for  the  mean  square  vertical  velocity 
becomes 


l 

J 


t  =  -/4 


cos 


2  2  TT  t 


t  =  0 


-2 

w 


2tt2  h2 


(22) 


The  equation  for  root  mean  .square  vertical  velocity  is  then 


(23) 


Expressions  for  h  in  terms  of  the  slope  of  an  isentropic  surface  and  for  in 
terms  of  stability  and  wind  shear  will  now  be  derived. 

In  a  turbulent  source  region,  energy  for  wave  motion  is  continually  being  produced 
at  a  constant  rats.  As  the  air  passes  downwind  beyond  this  region,  the  wave  energy 
may  be  gradually  dissipated  by  mixing  with  air  outside  the  turbulent  layer,  but  any 
such  dissipation  will  not  be  considered  here.  An  exoression  for  the  rate  of  pro¬ 
duction  of  wave  energy  per  unit  mass  can  be  defined  Dy 


_  t  =  7/4 

dE/dt  =  ^  J*  d/dt  (w2/2)  dt 
*t  =  0 
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Hence  the  resulting  wave  amplitude,  and  from  Eq.  (22)  the  mean  square  vertical 
velocity,  depend  upon  the  mean  rate  at  which  energy  becomes  available.  This  energy 
for  wave  motion  must  be  supplied  by  the  initial  perturbation. 

In  fluid  mechanics,  it  is  known  that  a  submerged  object,  such  as  a  weir,  may  dis¬ 
turb  the  flow  by  producing  a  standing  wave  over  the  object.  When  the  wave  ampli¬ 
tude  is  large  relative  to  the  depth  of  the  fluid,  the  familiar  hydraulic  jump  occurs. 
Upon  approaching  the  submerged  object,  the  flow  is  subjected  to  contraction  and 
acceleration.  This  is  followed  by  a  rapid  expansion  and  deceleration  in  the  standing 
wave  to  the  lee  of  the  obstruction,  e.g.,  see  King  (20),  Air  flowing  over  a  moun¬ 
tain  range  or  mound  of  cold  air  appears  to  behave  similarly,  e.g. ,  see  Queney 
et.  al.  (21). 

In  general,  any  change  In  conditions  of  flow  which  is  accompanied  by  contraction 
or  expansion  of  the  depth  of  fluid  produces  a  series  of  waves  on  the  sloping  fluid 
surface.  If  the  change  occurs  suddenly,  such  as  when  the  fluid  passes  over  a 
fixed  submerged  object,  an  abrupt  standing  wave  usually  occurs.  In  a  stably 
stratified  compressible  fluid  of  great  depth,  contraction  or  expansion  would  be 
associated  with  mass  divergence  and  hence  vertical  motion.  Assuming  adiabatic 
temperature  changes,  vertical  motion  would  produce  horizontal  temperature  gra¬ 
dients  and  hence  sloping  isentropic  surfaces.  Therefore,  energy  for  wave  forma¬ 
tion  might  be  a  function  of  the  mass  divergence  which  is  related  to  the  horizontal 
grrdients  of  both  density  and  velocity.  An  expression  which  takes  into  account  both 
mass  divergence  and  vertical  motion  is  the  pressure  tendency  equation,  Homboe, 
et  al.  (22)  which  may  be  written  as. 


(3P/3t)z  =  -g  j  Vh.(oV)6z  +  (gow)z 


(26) 
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In  Eq.  (26),  the  local  pressure  tendency  at  any  height  z  can  be  represented 
explicitely  as  the  change  in  weight  of  the  vertical  air  column  above.  This  change 
is  caused  in  part  by  the  horizontal  divergence  in  mass  above  the  level  and  in  part 
by  the  vertical  transport  of  mass  through  the  level.  The  horizontal  mass  diver¬ 
gence  is  defined  by 


7  .(pV)  =  a/dx(ou)  +  a/by(pv) 


(27) 


and  so  is  dependent  upon  the  horizontal  gradients  of  velocity  and  density. 

In  a  stationary  system,  such  as  a  mountain  wave,  the  horizontal  density  and  velocity 
gradients  are  fixed  and  the  local  pressure  tendency  term  in  Eq.  (26)  is  zero.  An 
expression  for  the  rate  of  energy  production  per  unit  mass  dE/dt  can  be  obtained 
by  dividing  the  remaining  two  terms  in  the  pressure  tendency  equation  by  the  mean 
density,  hence 


2 

VH*  (pV)  6z  =  gW 


z 


1 


(28) 


where  the  Integration  is  taken  over  the  thickness  of  the  layer  with  horizontal  density 
gradients,  i.  e. ,  the  layer  undergoing  vertical  motion.  The  inference  is  made  that 
whenever  a  perturbation  in  the  flow  causes  a  vertical  displacement  of  a  stably 
stratified  layer,  energy  for  wave  motion  is  produced,  and  according  to  Eq,  (28) 
the  mean  rate  at  which  this  energy  becomes  available  is  directly  proportional  to 
the  vertical  velocity  associated  with  the  initial  perturbation. 

For  an  adiabatic  process,  an  estimated  value  for  vertical  velocity  can  be  obtained 
from  the  product  of  the  horizontal  wind  speed  and  the  vertical  slope  of  the  isen- 
tropic  surfaces  in  the  direction  of  the  wind.  Of  course,  for  air  to  move  vertically 
along  such  a  sloping  surface,  it  must  be  moving  relative  to  the  surface,  i.e. , 
moving  faster  or  slower  than  the  apparent  foreward  motion  of  the  isentropic  sur¬ 
faces.  In  a  stationary  system,  such  as  a  standing  wave,  the  isotherms  remain 
stationary  and  the  sloping  isentropic  surfaces  define  the  trajectori?  s  of  the  air 
passing  through  the  system.  If  s  is  the  vertical  slope  of  an  isentropic  surface 
and  ur  is  the  component  of  the  mean  wind  speed  in  the  direction  of  the  slope,  then 
the  absolute  value  of  the  vertical  velocity  can  be  represented  by 


(29) 
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In  the  special  situation  where  the  isotherms  on  a  constant  pressure  level  are 
oriented  at  right  angles  to  the  height  contours  and  are  stationary  in  time,  then 
ur  is  simply  equal  to  the  mean  wind  speed  u  , 

Upon  substituting  this  expression  for  w  in  Eq.  (28),  we  find, 


dE/dt  =  g  j  uf  s 


where  ir  and  s  are  mean  values  taken  over  the  layer  of  interest.  Upon  equating 
this  expression  to  the  rate  of  wave  energy  production  in  Eq.  (25),  we  find 


s|  =  ^JL. 
t3 


(31) 


and  upon  rearranging  terms,  the  expression  for  amplitude  h  is 


h 


2 


(32) 


2 

Substituting  for  h  from  Eq.  (22) ,  the  expression  for  mean  square  vertical  velocity 
becomes 


(33) 


Substituting  for  t  from  Eq.  (18),  we  find 


(34) 


This  expression  states  that  for  a  given  set  of  initial  conditions,  i.e.,  given  values 
for  ur  ,  s  and  u  ,  the  ’  can  square  vertical  velocity  associated  with  the  wave  can 
be  expressed  as  a  function  of  wavelength  only. 
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An  expression  for  wavelength  in  terms  of  local  stability  and  wind  shear  will  now 
be  postulated.  In  a  stable  atmosphere,  the  frequency  v  is  defined  by 


v2  =  f  (y- vd)  =  f  ae/az 


(35) 


where  0  is  the  potential  temperature.  Substituting  from  Eqs.  (17)  and  (18),  the 
expression  for  wavelength  is  found  to  be 


2ttu 

^ g/e  ae/dz 


(30) 


However,  this  equation  does  not  take  into  account  the  possible  effects  of  vertical 
wind  shea".  It  is  well  known  that  shear  acts  to  restrict  the  vertical  dimensions 
of  turbulent  eddies  even  in  the  absence  of  stability.  One  method  of  introducing 
the  effect  of  shear  into  the  model  is  to  modify  Eq.  (35)  to  read 


(37) 


where  B2 


& 

0  yi 


The  term 


ou 

SI 


is  the  vector  wind  shear. 


Upon  substituting  Eq.  (37)  for  X  in  Eq.  (34),  the  expression  for  root  mean  square 
vertici  elocity  becomes 


u 


(38) 


The  slope  of  the  isentropic  surfaces  can  be  calculated  from  three  dimensional  fields 
of  temperature  and  pressure,  or  by  analyzing  the  height  contours  of  an  isentropic 
surface.  The  difficulty  here  is  the  tendency  to  equally  space  the  height  contours 
between  points  of  observation  where  actually  they  may  be  crowded  somewhere  in 
between.  The  thermal  wind  equation  may  prove  helpful  in  deciding  whether  the 
mean  temperature  gradient  through  a  layer  is  concentrated  near  a  point  of  obser¬ 
vation.  The  thermal  wind  equation  may  be  written 
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3  P/3n  =  3u/az 


(39) 


where  f  is  the  coriolis  parameter  and  equals  2  0  sin  cd  or  about  l(f2 3  4  at  mid¬ 
latitudes.  The  n-direction  is  a  horizontal  vector  oriented  peipendicularly  to  the 
wind  shear  vector.  The  expression  for  the  slope  is 


ap  ,  ae 

an  '  3z 


(40) 


and  of  course.  Up  in  Eq.  (38)  is  the  component  of  the  mean  layer  wind  in  the 
n-direction.  Substituting  Eq.  (39)  in  this  expression  for  S  ,  we  get 


8 


=  Au/AP 

o 


(41) 


where  Au  is  the  magnitude  of  the  vector  wind  shear  and  A  6  is  the  difference  in 
potential  temperature  through  the  same  vertical  layer,  e.  g. ,  a  layer  between  stan¬ 
dard  wind  reporting  levels. 

When  computing  RMS  vertical  velocities,  Eq.  (38)  could  be  used.  However,  in 
order  to  determine  whether  the  computed  value  is  reasonable,  X  may  be  computed 
from  Eq.  (37)  and  it  should  have  a  value,  on  the  average,  cf  a  few  thousand  meters. 
Next  S  can  be  computed  from  Eq.  (41),  or  in  areas  between  data  stations  from 
an  isentropic  height  analysis,  S  should  have  an  average  value  of  the  order  of  10"^. 
Then  a  value  for  h  can  be  obtained  from  Eq.  (32)  upon  substituting  t  =  X/u. 

It  should  average  around  one  to  a  few  hundred  meters.  Finally,  the  RMS  vertical 
velocity  can  be  obtained  from  Eq.  (23). 

Examples  of  World  Wide  Application  of  Both  Models 

In  order  to  test  the  applicability  of  both  models,  it  was  decided  to  analyze  data  for 
a  sample  of  the  HICAT  turbulence  tests  which  met  the  following  criteria: 

1.  The  radiosonde  and  rawinsonde  data  were  concurrent  with  the  aircraft 
flights. 

2.  The  observed  turbulence  occurred  over  or  very  near  a  radiosonde 
station. 

3.  The  tests  include  at  least  three  geographical  areas. 

4.  The  tests  include  observed  turbulence  intensities  ranging  from  light  to 
severe. 
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The  tests  chosen  are  listed  in  Table  m  along  with  the  particular  runs  and  geograph¬ 
ical  areas.  The  areas  are  Western  United  States,  Hawaii  and  Eastern  Australia. 

In  the  Australian  Test  102,  there  were  two  altitude  bands  with  turbulence  -  at 
60  K  ft  near  Brisbane  and  at  66  K  ft  near  William  at  own  -  so  both  are  treated  as 
separate  samples.  As  shown  in  Figures  29  and  30,  two  isentiopic  height  analyses 
were  made  corresponding  approximately  to  these  two  altitudes.  In  the  other  geo¬ 
graphical  areas,  the  slope  term  was  obtained  from  the  thermal  wind  equation. 

For  several  of  these  tests,  no  RMS  vertical  velocity  data  were  obtained.  However, 
RMS  incremental  eg  normal  acceleration  data  are  available.  These  data  are  based 
on  acceleration  of  the  aircraft  and  hence  do  not  include  the  component  representing 
the  acceleration  of  the  air  relative  to  the  aircraft.  Nevertheless,  the  acceleration 
of  the  aircraft  should  be  highly  correlated  with  the  intensity  of  turbulence  so  that 
computed  RMS  accelerations  based  on  the  proposed  model  should  correlate  well 
with  the  aircraft  acceleration  data  if  the  models  are  applicable.  In  the  previous 
section,  an  equation  was  derived  for  RMS  vertical  acceleration  which  an  aircraft 
would  experience  if  it  responded  completely  to  all  changes  in  vertical  motion  of 
the  air.  As  the  aircraft  response  is  usually  less,  the  computed  RMS  acceleration 
should  tend  to  be  high.  The  derived  Eq.  {15)  expressed  in  terms  of  the  acceleration 
of  gravity  is 


RMS  Aajj/g 


2^~2~  tt2h  uv 


■gx‘ 


(42) 


where  v  is  the  aircraft  speea  . 

s 

As  shown  in  Table  HI,  the  wavelengths,  amplitudes  and  RMS  Aajj/g  values  are 
computed  using  both  models.  The  observed  RMS  Aa^/g  values  are  given  in  the 
last  columns  In  terms  of  the  average  value  for  all  runs  in  the  area  and  the  range 
of  average  values  experienced  over  all  the  individual  runs. 


^Actually  v8  is  the  speed  of  the  aircraft  relative  to  the  wave  pattern.  In  a  station¬ 
ary  wave,  v8  is  the  ground  speed,  i.e.,  the  speed  relative  to  the  air  passing 
through  the  wave.  Since  the  speed  of  the  aircraft  is  always  much  greater  than 
the  wind  speed  and  in  most  instances  any  possible  progression  of  the  wave  pattern 
will  be  unknown,  any  correction  for  the  value  of  v8  does  not  seem  justified  at 
this  time.  Nevertheless,  it  may  account  for  the  occasional  differences  in  intensity 
of  turbulence  noted  by  pilots  between  downwind  and  upwind  runs.  These  differences 
were  especially  prevalent  in  Australia  where  the  wind  speeds  were  generally  greater 
than  in  other  areas. 
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Figure  29.  Isentroplc  Height  Analysis  for  0  =  450*A,  Southeastern  Australia, 
27  Jul  1966.  2300Z  (Flight  102) 
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Figure  30.  isentropic  Height  Analysis  for  6  «  490®A,  Southeastern  Australia, 
27  Jul  1966,  2300Z  (Flight  102) 
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TABLE  HI 


COMPARISON  OF  COMPUTED  AND  OBSERVED  RMS  AaN/g  VALUES 


Test  No. 
and 

Area 

Model 

Wave 

Length 

(M) 

*  ...  .  RMS 

Amplitude  A  , 

(M)  iaN/8 

Average  Observed  RMS  AaN/g 
Average  Range  over 

of  Runs  all  Runs 

77 

I 

7248 

213 

.043 

.036 

.027-.  050 

Hawaii 

n* 

4355 

144 

.081 

(Light) 

102 

I 

4933 

183 

.082 

.050-. 059 

Runs  12-14 

.056 

William  st  own 

Australia 

n** 

4605 

135 

.069 

(Light  to  Moderate 

114 

i 

2776 

147 

.140 

.051-. 078 

l  as  Vegas 

.060 

Nevada 

n* 

3665 

129 

.068 

(Moderate) 

102 

i 

1882 

69 

.115 

.038-. 095 

Runs  2-9 

.063 

Brisbane 

Australia 

n** 

3044 

85 

.060 

(Moderate) 

267  I 

3638 

178 

.169 

.  034- . 162 

Albuquerque 

.083 

New  Mexico  II* 

3851 

235 

.138 

(Severe) 

*Slope  determined  from  geostropic  thermal  wind  shear 
**Slope  determined  from  isentropic  height  analysis 


The  general  agreement  between  corresponding  computed  quantitatives  for  the  two 
methods  is  good  considering  the  fact  that  the  two  models  involve  entirely  different 
data.  In  Model  I,  the  data  are  obtained  directly  from  the  radiosonde  temperature- 
height  curve;  while  in  Model  II,  the  computations  are  based  upon  an  assumed  rate 
of  vertical  ascent  along  a  sloping  isentropic  surface. 

In  general,  it  appears  that  the  range  of  values  calculated  by  Model  I  yields  ItMS  Aa^/g 
values  which  correlated  slightly  better  with  the  indicated  intensities  of  turbulence, 
i.e. ,  light  to  severe.  They  are  also  greater  than  the  average,  as  was  expected, 
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but  within  a  factor  of  two  of  the  observed  value  for  the  most  turbulence  run  in  each 
area.  On  the  other  hand.  Model  n  yields  RMS  A  a^/g  values  which  are  in  general 
closer  to  the  observed  average  values  of  all  the  runs  in  each  area  but  generally 
less  than  the  observed  maximum  average  value. 

STATISTICAL  STUDIES 

The  literature  describing  the  results  of  statistical  and  case  history  studies  of  rela¬ 
tionships  between  meteorological  conditions  and  high  altitude  clear  air  turbulence 
is  summarised  in  Table  I,  pg„  3-5.  Generally,  the  meteorological  conditions  have  been 
expressed  explicitly  as  functions  of  the  atmospheric  temperature,  pressure,  and 
wind  fields.  The  temperatures,  pressures,  and  winds  at  HICAT  flight  levels  are 
measured  at  stations  that  are  spaced  two  nr  more  hundreds  of  miles  apart.  Temper¬ 
atures,  heights  of  pressure  surfaces,  and  winds  are  reported  for  specific  intervals. 

The  observations  are  made  at  00Z  and  12 Z  time  except  in  some  areas  in  which  only 
the  00Z  observations  are  made.  Horizontal  and  vertical  temperature  gradients 
ami  wind  shears  are  derived  quantities.  With  few  exceptions  the  horizontal 
temperature  gradients  and  the  horizontal  wind  shears  are  obtained  by  drawing 
isolines  on  charts.  This  method  introduces  a  subjective  element  into  the  analysis. 

In  the  following  discussion  an  attempt  has  been  made  to  keep  the  subjective  element 
to  a  minimum. 

The  criteria  used  for  selecting  the  cases  for  study  are  as  follows: 

A  "run"  as  given  by  Crooks,  floblit  and  Prophet  (2)  was  designated 
as  a  sample  of  turbulence  if  the  run  was  within  100  miles  of  a  RAOJB 
or  RAWIN  station.  A  sample  was  designated  "non-turbuient"  if  no 
"’run"  was  shown  within  150  miles  of  a  station  and  the  flight  track  was 
within  100  mile*,  of  a  station. 

All  such  'turbulence"  and  "no  turbulence"  cases  from  HICAT  flights  54-179  were 
used  in  the  study.  These  flights  were  made  in  Australia,  New  Zealand,  Hawaii, 

Alaska,  California,  Massachusetts  and  Puerto  Rico.  Frequency  distributions  for 
eight  functions  of  the  temperature  and  wind  fields  are  shown  in  Figures  31  through 
32.  The  interval  limits  were  choser  to  produce  approximately  the  same  number 
of  cases  in  each  interval.  Chi-square  tests  with  four  degrees  of  freedom  and 
Kolmogorov-Smirnov  tests  were  used  to  indicate  the  degree  to  which  the  turbulence 
and  no  turbulence  cases  may  be  considered  to  come  from  differing  populations. 
x  ■  9.49  is  the  critical  value  for  x  with  5  degrees  of  freedom.  Degrees  of 
freedom  are  determined  by  the  number  of  class  intervals.  For  a  y2  value  cal¬ 
culated  from  the  observed  data  greater  than  9. 49,  the  hypothesis  that  the  two  fre¬ 
quencies  are  from  the  same  population  is  rejected  at  the  95%  confidence  level. 

In  the  Kolmogorov-Smirnov  test,  K-S  =  1. 36  is  the  critical  value  above  which 
the  null  hypothesis  that  the  two  frequencies  come  from  the  same  population  is 
rejected  at  the  5%  level.  If  the  critical  statistic  obtained  from  the  observed  data 
is  less  than  1. 36,  the  null  hypothesis  is  accepted  (the  two  frequencies  are  assumed 
to  come  from  the  same  distribution).  The  results  of  the  two  tests  differ  mainly 
because  of  the  greater  power  of  the  Kolmogorov  -Smirnov  test  which,  unlike  the 
Chi-square  test,  considers  sample  Bize,  The  tests  are  an  estimate  of  by  how  much 
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the  two  distributions  differ.  In  both  the  Chi-square  and  Kolmogorov -Smirnov  tests 
larger  numbers  for  the  critical  value  indicate  less  chance  that  the  distributions 
were  drawn  from  the  same  population. 

The  results  obtained  by  applying  the  Chi-square  and  die  Kolmogorov-Sm irnov  tests 
are  presented  in  Table  IV. 


TABLE  IV 

CHI-SQUARE  (X2)  AND  KO LMOGOROV-SMIRNOV  (K-S)  TESTS  APPLIED  TO 
METEOROLOGICAL  VARIABLES  FOR  HICAT  FLIGHTS  54  THROUGH  179 


x2 

(Deg.  Freedom:  4) 

K-S 

Critical  Values  p™1»b^  - 

Probability  *  .99 

9.49 

13.30 

1.36 

1.63 

V 

Mean  scalar  wind  speed  (5000  to 

10, 000  ft  layers) 

28. 8 

1.08 

AV/Az 

Vertical  scalar  wind  shear 

86.3 

2.30 

Iav|/Az 

Magnitude  of  vertical  vector  wind 
shear 

126.5 

2.69 

V(lA?|/Az) 

Product  mean  scalar  wind  speed 
and  vertical  vector  wind  shear 

89.8 

2.27 

5T/dz 

Vertical  temperature  gradient  at 
flight  level 

76.8 

1.94 

(AT|/Az 

Absolute  value  of  temperature  change  65. 1 
(layer  +  2000  ft  of  flight  level) 

2.05 

( oT/6z)* 

Min  vertical  temperature  gradient 
(layer  +  2000  ft  of  flight  level) 

90.4 

2.11 

Ri 

f  cae/az)(lAv!/Azf2,  ae/azis 
min  +  2000  ft  of  flight  level 

144.1 

2.89 

Figure  31a  and  the  first  rcw  in  Table  IV  indicate  that  at  the  altitudes  of  the  HICAT 
flights  no  significant  distinction  can  be  made  between  the  population  of  turbulence 
and  non-turbulence  cases  for  scalar  winds.  The  scalar  winds  at  the  flight  altitudes 
were  obtained  by  interpolating  between  the  winds  reported  by  the  meteorological 
station.  This  was  necessary  because  the  altitudes  for  reporting  the  winds  do  not 
correspond  with  the  flight  altitudes.  The  data  presented  in  Figure  31a  does. 
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iMX?arer;  justify  further  study  with  a  ’arger  sample  for  the  strong  wind  cases. 

See*  e.g.,  Briggs  and  Roach  (23). 

Figures  31b  and  c  and  rows  2  and  3  of  Table  IV  indicate  the  vertical  vector  or  the 
scalar  wind  shear  may  be  used  to  separate,  to  a  significant  degree,  the  populations 
of  turbulence  and  non-turbulence.  The  use  of  the  vector  wind  shear  produces  the 
greater  success  in  the  separation  of  the  populations.  This  conclusion  is  derived 
from  the  values  in  Table  IV  for  the  Chi-square  test  and  the  Komolgorw -Smirnov 
teat.  In  SI  cases  of  the  vertical  vector  shear  equal  to  or  less  than  2  kts/1000  ft 
7  were  turbulent  and  44  non-turbulent.  There  were  twenty  instances  of  the  vertical 
vector  shear  equal  to  or  greater  than  5  kts/1000  ft.  Of  these,  eighteen  were  cases 
of  turbulence. 

The  product  of  the  scalar  wind  and  the  vertical  vector  shear  (Figure  3 Id  and  row 
4  of  Table  IV\  has  leas  value  as  an  aid  for  the  separation  of  turbulence  and  non¬ 
turbulence  populations  than  the  vector  wind  shear  alone. 

Several  methods  were  employed  in  analyzing  the  vertical  temperature  gradient  rela¬ 
tionship  to  turbulence: 

a)  Toe  vertical  gradient  at  turbulence  level  (Figure  32a); 

b)  The  absolute  value  of  the  average  gradient  over  a  layer  +  2000  ft  above 
and  below  the  level  of  turbulence  (Figure  32b); 

c)  The  minimum  gradient  in  any  layer  within  +  2000  ft  of  the  level  of  turbu¬ 
lence  (Figure  32c). 

Results  suggest  a  relationship  exists  between  large  decreases  in  actual  temperature 
f small  potential  temperature  gradients)  with  height  and  turbulence  (Figure  32a), 
although  several  cases  of  turbulence  were  accompanied  by  large  temperature 
increases  (stable  layers).  Only  3%  erf  the  non-turbulent  cases  had  vertical  temper¬ 
ature  gradients  <  -I’C/IOOO  ft,  compared  with  L9%  in  turbulence.  Changes  greater 
than  +1*C/1000  ft  were  associated  with  14%  of  the  non-tu^bulent  and  22%  of  the  tur¬ 
bulent  cases.  There  is  a  definite  trend  towards  die  association  of  large  fluctuations 
in  vertical  temperature  gradients  with  turbulence.  Conversely,  uniform  temperature 
gradients  are  more  likely  to  be  accompanied  by  smooth  flying  conditions. 

Richardson's  Number  (Figure  32d)  may  be  used  relatively  successfully  to  distinguish 
between  turbulence  and  non-turbulence.  Turbulence  was  recorded  in  32  of  42  cases 
with  Ri  <  15.  It  is  apparent,  however,  that  small  values  of  Ri,  recommended  by 
various  authors  as  the  necessary  criterion  for  the  occurrence  erf  turbulence  in  the 
troposphere,  are  rarely  observed  in  the  stratosphere.  There  was  but  1  case  in 
115  of  Ri  <  1,  a  turbu’ence  criterion  used  by  Briggs  and  Roach  (23)  and  Endlich 
and  Mancuso  (24)  but  rejected  by  Reiter  and  Lester  (25)  due  to  lack  to  lack  in 
quality  of  wind  and  temperature  measurements  in  the  stratosphere. 
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Figure  31.  Frequency  Distribution  of  Meteorological  Variables  for  Turbulent 
and  Nonturbulent  Cases 


PERCENTAGE  PERCENTAGE 


THE  DEVELOPMENT  OF  A  PRACTICAL  PROCEDURE  FOR  FORECASTING 
HIGH  ALTITUDE  CLEAR  AIR  TURBULENCE 

The  problem  of  developing  practical  methods  of  forecasting  high  altitude  clear  air 
turbulence  has  been  approached  from  three  points  of  view;  theoretical,  statistical, 
and  that  of  a  synoptic  meteorologist.  The  theoretical  work  was  presented  in 
Section  V  and  the  statistical  work  in  Section  V  .  In  March  1967  a  meteorologist 
was  assigned  to  the  WU-2  flight  crew.  This  assignment  provided  an  excellent 
opportunity  for  a  meteorologist  to  work  closely  with  the  pilots  in  selecting  the 
HICAT  search  areas,  the  search  procedures,  and  to  modify  the  search  plan  in 
mid-flight  whenever  this  appeared  to  be  desirable.  In  addition,  the  meteorologist 
obtained  direct  practical  experience  in  operating  under  field  conditions. 

The  method  of  preparing  the  forecasts  of  the  HICAT  have  changed  as  the  program 
progressed.  These  changes  have  been  associated  with  differing  communication 
and  meteorological  organizations  at  the  various  bases  of  operation  and  with  meteoro¬ 
logical  conditions  and  added  experience.  The  present  state  of  development  o'  a 
forecasting  procedure  is  illustrated  by  two  case  histories  described  In  detail 
below. 

Case  of  Severe  Clear  Air  Turbulence 


Synoptic  Weather  Patterns 

The  weather  pattern  generally  believed  to  favor  the  production  of  CAT  below  the 
tropopause  is  typically  that  in  which  a  well  established  low  pressure  center  near 
the  500  mb  level  moves  across  the  western  states.  The  polar  jet  stream  dips 
far  to  the  south,  bringing  cold  air  at  most  levels  from  the  surface  to  the  tropopause. 
If  the  trough  persists  through  the  tropopause,  the  temperatures  in  the  stratosphere 
over  the  southern  station  portion  increase.  At  the  70  n  ,  iovel  alternate  north- 
south  zones  of  relatively  hign  and  low  temperatures  are  observed. 

Synoptic  patterns  such  as  these,  which  are  primarily  a  winter  and  spring  pheno¬ 
mena,  seem  to  be  related  to  conditions  that  produce  abundant  clear  air  turbulence 
in  both  the  troposphere  and  stratosphere , 

These  typically  CAT  producing  synoptic  patterns  were  observed  to  occur  during 
two  different  periods  in  November  and  December  over  the  western  United  States. 

The  first  period  was  from  12  to  23  November  1967  rnd  the  second  from  29  November 
to  2  December  1967. 

Moderate  to  severe  CAT  was  encountered  at  approximately  50, 000  feet  over  Santa 
Barbara,  California,  during  flight  number  257  on  17  November  1967.  A  post  flight 
analysis  of  this  mission  provided  clues  that  were  utilized  in  refining  forecast  pro¬ 
cedures.  These  procedures  were  applied  at  the  beginning  of  the  next  favorable 
turbulent  period.  Results  of  the  HICAT  flights  during  this,  period  were  as  follows: 
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Test 

No. 

Date 

Place 

Turbulence 

Rating 

Incremental  c.g. 
Acceleration 

262 

29  Nov  1967 

Albuquerque,  N.M. 

Moderate 

(±.40g) 

264 

30  Nov  1967 

Grand  Junction,  Colo. 

Severe 

(+.65g) 

265 

30  Nov  1967 

Grand  Junction,  Colo. 

Severe 

(+.65g) 

(night) 

266 

1  Dec  1967 

Albuquerque,  N.M. 

Severe 

(+  1.  Og,  -0. 7g) 

267 

1  Dec  1967 

Albuquerque,  N.M. 

Severe 

(i.60g) 

Forecast  Procedure 

Prior  to  the  beginning  of  HICAT  operations  on  16  November  1967  from  Edwards 
Air  Force  Base,  California,  a  number  of  meteorological  analyses  techniques  had 
been  tested  with  varying  degrees  of  success.  Analyses  of  tropopause  slopes,  hori¬ 
zontal  and  vertical  wind  shear,  Richardson  Numbers,  isobaric  surface,  pressure 
height  changes,  and  thickness  patterns  did  not  indicate  a  pattern  of  consistency 
insofar  as  a  reliable  forecast  technique  was  concerned.  One  notable  feature,  how¬ 
ever,  did  stand  out  in  its  relation  to  turbulence  and  this  was  a  distinctive  isothermal 
pattern.  With  this  idea  in  mind  it  was  decided  to  analyze  three  surfaces  (100,  70, 
and  50  mb)  for  isotherms  and  wind  streamlines  only.  These  charts  indicated  that 
the  gradient  of  temperature  along  the  constant  pressure  surface  varied  to  a  signifi¬ 
cant  degree. 

It  can  be  stated  with  a  reasonable  degree  of  confidence,  that  CAT  will  occur  in  an 
area  where  the  temperature  gradient  is  large.  The  critical  value  of  the  hoiizontal 
temperature  gradient  appears  to  be  about  5°C  per  120  rnn.  It  is  significant  that 
the  Meteorological  Department  of  Trans  World  Airlines  in  their  Technical  Bulletin 
Number  61-2,  observed  such  a  horizontal  temperature  gradient  to  be  associated 
with  CAT  in  the  tropospheric  range  of  15, 000  -  40,  000  feet.  For  gradients  less 
than  this  turbulence  should  be  light  to  non-existent.  For  gradients  higher  than 
this,  moderate  or  severe  CAT  may  be  expected. 

This  is  considered  a  good  generalization  but  should  be  qualified  to  some  extent. 

For  example,  it  cannot  be  stated  unequivocally  that  these  criteria  would  apply  to 
CAT  associated  with  convective  activity  or  a  mountain  wave  although  there  may 
very  well  be  a  similar  relation. 

Constant  Pressure  Surface  Analysis 

Figure  33  shows  the  route  of  the  flight  made  on  1  Dec  1967.  The  pilot  flew  to 
Albuquerque,  New  Mexico,  located  the  turbulence  directly  over  Albuquerque  and 
began  a  vertical  search  pattern  that  a  short  time  later  resulted  in  the  location  of 
severe  CAT.  Figure  34  is  the  pattern  at  the  70  mb  level  that  was  associated  with 
the  severe  turbulence. 
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The  analysis  of  the  70  mb  level  in  Figure  34  shows  the  strong  horizontal  temper¬ 
ature  gradient  that  exists  between  Albuquerque,  New  Mexico,  and  El  Faso,  Texas, 
at  00Z  2  December  1967.  These  RAOB's  were  being  taken  at  the  same  time  severe 
CAT  was  encountered  over  Albuquerque  on  Test  266.  The  horizontal  temperature 
gradient  on  the  constant  pressure  surface  of  70  mb  was  approximately  10#/180  nm. 
This  is  considerably  above  the  critical  value  of  5 ®/120  nm  mentioned  above.  It 
should  be  noted  that  the  70  mb  chart  in  Figure  34  was  not  that  used  for  the  forecast 
for  Test  266  but  was  the  analysis  nearest  the  actual  time  the  aircraft  was  over 
Albuquerque.  Twelve  and  twenty-four  hour  analyses,  prior  to  the  00Z  chart  did 
show,  however,  a  steady  progression  and  intensification  of  the  horizontal  temper¬ 
ature  gradient  across  Utah,  Arizona  and  New  Mexico. 


Figure  35  shows  a  plot  of  temperature  as  a  function  of  pressure  over  Grand  Junction, 
Colo. ,  Albuquerque,  New  Mexico,  and  El  Paso,  Texas  for  00Z  2  Dec  1967  and 
corresponding  to  the  time  of  the  horizontal  70  mb  surface  analysis  (Figure  34). 

The  irregular  appearance  of  the  sounding  is  quite  obvious  as  compared  to  that  of 
Grand  Junction,  Colo. ,  which  Figure  34  shows  to  be  in  an  area  of  small  horizontal 
temperature  gradient.  The  turbulence  over  Albuquerque  was  reported  by  the  pilot 
to  be  light  to  severe  at  altitudes  51, 000  to  57, 000  feet.  The  altitudes  reported 
were  those  read  directly  from  the  altimeter  which  was  set  at  29. 92  inches. 


The  fact  that  CAT  is  associated  with  a  strong  temperature  inversion  been 
recognized  for  some  time.  S.  M.  Serebreny  (26)  stated  "Turbulence  is  a  result 
of  an  abrupt  increase  or  decrease  of  wind  velocity  in  the  thermal  gradient.  The 
stronger  the  vertical  shear  the  more  t  table  the  lapse  rate".  Frederick  B.  Tlaymond 
( 19) »  in  evaluating  953  forecasts  for  110  U-2  flights  out  of  Davis  Monthan  Air 
Force  Base,  Arizona,  observed  thaf  excellent  forecast  verification  was  achieved 
by  predicting  CAT  to  occur  in  areas  and  near  levels  of  inversions  as  follows: 


1. 5°C/1000  ft 
1.5°C  to  2„5°C/100G  ft 
2.5°C  to  4. 0°C/1000  ft 
4. 0°C/1000  ft 


No  CAT 
Light  CAT 
Moderate  CAT 
Severe  CAT 


Isentropic  Analysis 

Potential  temperature  is  a  function  of  pressure  and  temperature,  hence  charts  of 
surfaces  of  equal  potential  temperature  (isentropic  surfaces)  are  simply  an  alter¬ 
nate  method  of  indicating  temperature  gradients.  Often  isentropic  cross-sections 
are  the  most  convenient  method  of  determining  the  extent  of  baroclinicity.  Reiter 
and  Hayman  (27) ,  commenting  in  the  paper  "The  Nature  of  Clear  Air  Turbulence", 
stated  'It  became  evident,  however,  that  all  cases  of  moderate  and  severe  CAT 
were  located  in,  or  very  close  to,  the  axis  of  a  drop  of  the  isentropic  surfaces  - 
the  isentrope  trough". 

Figure  36  is  an  isentropic  cross-section  through  Salt  Lake  City,  Utah,  Grand 
Junction,  Colorado,  Albuquerque,  New  Mexico,  and  El  Paso,  Texas,  as  of  00Z 
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1SENTR0PIC  SURFACE  00 Z 


Figure  36.  Isentropic  Analysis,  2  De  s  1967,  0000Z 


2  Dec  1967.  It  corresponds  to  the  flight  track  in  Figure  33  and  the  analyses  in 
Figures  34  and  35. 


The  baroclinic  zones  are  quite  pronounced,  particularly  over  Albuquerque  and 
El  Paso.  As  stated  above  by  Reiter  and  Hayman,  the  turbulence  was  encountered 
near  the  axis  of  an  i  sen  trope  trough.  It  should  be  noted  however,  that  the  turbu¬ 
lence  had  decreased  to  very  light  at  57, 500  feet  which  was  at  the  low  point  of  the 
isentrope  trough.  The  severe  CAT  was  observed  to  be  at  about  54, 000  feet  which 
was  slightly  below  the  level  of  the  inversion.  Although  the  turbulence  was, in  fact, 
along  the  axis  of  the  isentrope  trough  the  severe  turbulence  was  actually  in  an 
isentropic  hump.  It  is  perhaps  more  significant  that  the  turbulence  occurred  near 
the  "ertical  axis  of  a  zone  in  which  there  was  considerable  packing  and  spreading 
of  the  isentrope  s. 

One  of  the  most  difficult  parameters  to  forecast  in  regard  to  clear  air  turbulence 
is  the  altitude  at  which  it  is  likely  to  occur.  If  the  turbulent  area  moves  along, 
or  close  to  the  isentropic  surface  then  the  isentropic  analysis  should  prove  to  be 
an  aid  in  forecasting  levels  of  turbulence. 

No  Turbulence  Case 


It  is  interesting  to  note  that  the  bibliography  in  this  report  lists  many  titles  of 
studies  concerning  clear  air  turbulence  but  no  titles  of  studies  devoted  to  smooth 
air  or  lack  of  turbulence. 

To  be  able  to  closely  define  the  limits  of  any  problem  is,  of  course,  a  distinct 
advantage.  With  this  thought  in  mind  numerous  cases  where  CAT  was  not  found 
have  been  analyzed.  It  can  be  stated  that  there  are  certain  characteristics  of 
synoptic  patterns  that  are  associated  with  conditions  of  no  turbulence. 

Figure  37  shows  the  route  of  flight  on  30  June  1967.  From  Figure  38  it  can  be  seen 
that  the  flight  at  the  70  mb  level  was  in  an  area  of  weak  horizontal  temperature  gra¬ 
dient.  The  -50*C  and  -52 "C  isotherms  are  approximately  300  nm  apart.  The  iso¬ 
therms  are  nearly  parallel  to  the  streamlines. 

Figure  39  shows  the  RAOBS  at  Moossonee,  Ont. ,  Maniwaki,  Que. ,  and  Caribou, 
Maine.  The  soundings  are  relatively  straight  as  compared  to  those  in  Figure  35 
indicating  small  temperature  changes  in  the  vertical. 

The  isentropic  analysis  in  Figure  40  shows  a  smooth  isentropic  surface  with  a  alight 
downward  tilt  towards  the  colder  air  to  the  north. 

The  obviously  smooth  surfaces  of  30  June  1967  may  be  compared  to  the  relatively 
steeply  sloping  ioeritropic  surface  of  2  Dec  1967.  In  short  the  analyses  of  30  June 
1967  reflect  a  pattern  of  weak  thermal  activity.  Enough  cases  v.*ve  been  examined 
to  indicate  that  when  upper  level  synoptic  conditions  are  similar  to  those  in 
Figures  37,  38,  39,  and  40,  clear  air  turbulence  w:>  :  -  very  light  to  nonexistent. 

In  summary  then  it  may  be  said  that,  of  the  108  cases  examined  since  March  1967, 
if  the  flight  of  30  June  1967  is  typical  of  that  of  no-CAT  and  the  flight  of  1  Dec  1967 
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Figure  40.  Isenlropio  Analysis,  30  Jun  1967,  1200Z 


is  representative  of  conditions  associated  with  severe  to  extreme  clear  air  turbu¬ 
lence  significant  progress  has  been  made  towards  defining  the  limits  of  the  problem 
of  identifying  meteorological  parameters  associated  with  turbulence. 

HORIZONTAL  TEMPERATURE  GRADIENTS  IN  REGIONS  OF  CLEAR  AIR 
TURBULENCE 

Kadlec  (28,  29),  George  (30  ),  and  Helvey  (31 )  have  discussed  measure¬ 
ments  of  temperature  changes  that  occurred  in  regions  of  clear  air  turbulence. 

Most  of  these  measurements  were  made  at  the  flight  altitudes  of  commercial  jet 
aircraft  (25, 000  to  35, 000  ft).  A  few  of  the  measurements  were  made  in  the  50, 000 
to  60, 000  ft  range  of  altitude.  These  measurements  indicated  that  significant  clear 
air  turbulence  is  strongly  correlated  with  temperature  changes  that  often  are  as 
large  as  3TJ  in  2  minutes  of  flight  time. 

Crooks,  Hob  lit,  and  Prophet  (2)  presented  approximately  930  minutes  of  time 
histories  of  ambient  temperatures  that  are  synchronous  with  the  turbulence  observed 
in  HICAT  flights  No.  54  to  No.  178.  An  additional  543  minutes  of  time  histories 
were  available  for  the  present  study.  The  time  histories  of  the  ambient  temperatures 
were  displayed  by  plotting  12. 5  points/sec.  The  sensitivity  of  the  temperature 
measuring  device  is  such  that  it  does  not  adequately  follow  changes  of  temperature 
whose  frequencies  are  greater  than  one  cycle/sec.  In  the  study  described  in  the 
following  paragraphs  frequencies  of  temperature  changes  greater  than  one  cycle 
per  five  seconds  were  eliminated  by  smoothing. 

The  study  of  the  atmosphere  temperature  changes  associated  with  clear  air  turbu¬ 
lence  may  be  divided  conveniently  into  two  parts:  The  temperature  changes  observed 
within  a  region  cf  turbulence,  and  the  temperature  changes  observed  within  a  few 
minutes  preceding  the  entry  into  the  turbulent  region.  This  section  will  consist 
largely  of  a  discussion  of  the  temperature  changes  within  the  turbulent  region  only. 
Few  data  relevant  to  the  non-turbulent  region  were  received  in  time  to  be  included 
in  this  report. 

First,  the  definition  of  a  turbulent  region  given  by  Crooks,  Hoblit,  and  Prophet 
(2)  was  accented.  Their  designation  by  flight  number  and  run  number  was  also 
accepted.  Small  sample  calculations  of  the  correlation  coefficient  between  the 
temperature  changes  and  the  RMS  values  of  the  gust  velocities  (Figure  41a  and  b) 
indicated  that  the  correlation  coefficient  was  low.  Further,  the  changes  in  air¬ 
craft  altitude  during  the  •’run"  made  it  difficult  to  distinguish  between  horizontal 
and  vertical  temperature  gradients. 

Selected  time  histories  that  illustrate  the  relationship  between  clear  air  turbulence 
and  gust  velocities  are  shown  in  Figures  42  to  46.  A  comparison  of  two  time  his¬ 
tories  obtained  from  HICAT  flights  over  the  Cobar,  Australia  urea  is  shown  in 
Figure  42a  and  b.  One  without  turbulence  (flight  100,  run  9)  shows  nearly  isothermal 
conditions  and  one  with  moderate  turbulence  (flight  100,  run  6)  shows  short  time 
interval  temperature  changes  of  2.5°C  accompanying  the  turbulence.  The  runs 
are  separated  by  a  time  interval  of  1  1/2  hours. 
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Histories  of  Gust  Velocity  and  Temperature  Variations  for  Flight  100,  Runs  6  and 


Histories  of  Gust  Velocity  and  Temperature  Variations  for 
;  202,  Run  7  and  Flight  102,  Run  9 


The  variability  of  temperature  preceding  turbulence  contrasted  to  that  within  the 
turbulence  is  illustrated  for  flight  107,  run  8,  by  Figure  43.  These  changes 
occurred  above  relatively  flat  terrain  between  Laverton  and  Cobar  in  southeast 
Australia.  The  temperature  scale  has  been  enlarged  to  emphasize  changes  during 

turbulence. 

The  largest  temperature  change  for  any  flight  occurred  during  flight  265,  run  37 
(Figure  43b)  when  an  11.5*C  change  in  4  miles  was  recorded  along  with  severe 
turbulence.  The  flight  was  over  mountainous  terrain  near  Albuquerque,  New  Mexico. 

Although  statistical  studies  have  suggested  a  strong  relationship  between  turbulence 
and  atmospheric  temperature  changes,  the  characteristic  of  the  change  varies 
considerably  for  any  given  case.  For  example.  Figure  44a  shows  a  temperature 
change  of  9“C  in  3  miles  recorded  during  flight  202,  run  7  above  thunderstorm 
activity  northeast  of  Little  Rock,  Arkansas.  Flight  102,  run  9(Figure  44b)  shows 
a  gradual  rise  in  temperature  over  a  distance  of  several  miles,  although  a  drop 
in  the  temperature  erf  5T3  in  2  1/2  miles  appears  near  the  end  of  the  run.  The 
moderate  turbulence  accompanying  this  flight  was  at  61, 000  ft  on  the  lee  side  erf 
the  Great  Dividing  Range  in  southeast  Australia.  There  was  a  strong  jet  at  36, 000ft 
with  winds  exceeding  150  knots.  This  case  illustrates  that  the  effectiveness  of  a 
remote  sensing  devcie  may  depend  upon  the  relative  flight  direction.  Had  the 
aircraft  flown  at  180°  to  the  actual  heading  the  sharp  temperature  change  would 
have  probably  preceded  the  turbulence.  It  is  /■easonable  to  assume  this  type  of 
change  could  be  detected  more  readily  than  a  gradual  one. 

The  temperature  variations  for  flight  114,  run  13  and  run  14  (Figure  45a  and  b) 
east  of  the  ridge  line  of  the  Sierra  Nevada  Mountains  in  California,  are  indicative 
of  a  mountain  wave  situation.  The  heading  during  flight  114,  run  13  was  towards 
the  ridge  and  away  from  the  ridge  for  flight  114,  run  14.  Both  runs  were  at 
right  angles  to  the  ridge,  nearly  overlapping,  and  separated  by  only  a  few  minutes. 
Time  histories  for  flight  114,  run  13,  are  plotted  in  reverse  to  correspond  with 
the  wind  flow  which  is  west-southwest  and  perpendicular  to  the  ridge  both  at  moun¬ 
tain  top  and  flight  levels. 

A  trough- ridge  structure  in  the  horizontal  temperature  field  is  discernible  for  both 
runs.  The  distance  between  crests  is  approximately  9  to  10  miles,  with  the 
exception  of  the  first  wave  which  has  a  wavelength  about  three  times  as  large. 

In  a  detailed  report  on  1964  U-2  measurements  in  the  same  general  area,  Helvey 
(1967)  concluded  that  the  temperature  troughs,  or  cold  zones,  were  associated 
with  crests  in  the  mountain  wave.  Turbulence  was  found  to  be  concentrated  in 
the  zone  between  the  warm  axis  and  the  cold  axis,  the  warm  axis  being  oriented 
upwind  from  the  cold.  The  present  HICAT  flights  do  not  supply  any  conclusive 
evidence  to  verify  this.  There  is,  however,  a  case  of  severe  turbulence  in  the 
Denver,  Colorado  area,  flight  280,  run  10  (Figure  46)  where  the  greatest  gusts 
occurred  in  the  zone  designated  by  Helvey.  Figure  265,  run  37  (Figure  43b) 
also  verifies  this. 

The  examples  discussed  in  the  preceding  paragraphs  illustrate  that  the  gust  velocities 
are  not  simple  functions  of  the  variations  in  the  ambient  temperature.  The  flights 
through  the  turbulent  regions  designated  by  Crooks,  Hoblit  and  Prophet  (2)  had 
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time  durations  that  varied  from  20  seconds  to  15  minutes.  In  the  present  study 
only  those  turbulence  cases  that  had  a  duration  of  one  minute  or  more  were  studied. 
Various  methods  of  determining  the  degree  of  correlation  between  any  or  all  of 
the  components  of  the  gust  velocities  and  the  variation  of  temperature  were  tried. 
The  method  that  yielded  the  highest  correlation  coefficient  was  as  follows: 

a)  Each  run  was  used  as  one  sample.  The  largest  temperature  change 
for  any  20  second  period  during  which  the  altitude  did  not  change  by 
more  than  100  ft  was  determ  ned. 

b)  The  largest  change  in  gust  velocity  during  any  10  second  period  for 
each  of  the  three  components  was  determined. 

c)  The  temper  ature  change  and  the  gust  velocity  changes  need  not  be 
synchronous  nor  have  a  fixed  time  lag. 

2  2  2  1/2 

d)  The  quantity  (AUp  +  Au£  +  AUy  ^xnax  was  computed. 

e)  The  correlation  coefficient  between  this  quantity  and  the  temperature 
change  described  above  was  computed. 

The  results  of  the  computations  are  illustrated  in  Figures  47a  and  b.  For  73  casee 
selected  from  16  flights  from  the  records  from  flights  54  to  178  the  correlation 
coefficient  was  0. 76 .  For  69  cases  from  21  flights  from  flights  180  to  257  the 
correlation  coefficient  was  0. 82.  In  only  10  of  52  cases  was  the  maximum  temper¬ 
ature  change  less  than  1.5°C  if  (AU~  +  AU^  +  Aui;)  *'2  was  greater  than  or 

i;  Xj  V  llloA 

equal  to  25  ft/scc,  whereas  78  of  88  cases  had  temperature  changes  less  than 
l»5°Cwhen  (AU2  +  AU^  +  AUy  2  was  less  than  25  ft/sec. 

GEOGRAPHIC  DISTRIBUTION  OF  TURBULENCE 

In  previous  sections,  evidence  is  presented  that  high  altitude  clear  air  turbulence 
is  a  function  of,  in  part,  relatively  large  horizontal  temperature  gradients  at  the 
flight  level.  The  direct  determination  of  the  frequency  of  occurrence  of  these 
large  horizontal  temperature  gradients  for  all  the  principal  regions  of  the  earth 
was  beyond  the  scope  of  work  for  this  report.  However,  an  alternative  approach 
to  the  problem  was  feasible.  This  approach  was  based  upon  the  hypothesis  that 
regions  of  strong  horizontal  temperature  gradients  are  also  regions  in  which  the 
24  hour  temperature  changes  are  relatively  large.  This  presented  ah  opportunity 
to  use  the  Northern  Hemisphere  upper  atmosphere  temperature  data  collected 
during  the  International  Years  of  the  Quiet  Sun  (1964-1965).  These  data  were 
available  on  magnetic  tape.  Data  from  the  100  mb  level  for  the  western  portion 
of  the  Northern  Hemisphere  were  selected  for  analysis.  In  addition,  the  27  day 
period  of  V/U-2  HICAT  flight  operations  over  southeastern  Australia  were  included 
in  the  analysis.  The  100  mb  level  averages  close  to  53, 000  ft. 

For  purposes  of  the  analysis  the  Northern  Hemisphere  data  were  divided  by  season 
defined  as  winter  (Jan-Mar),  spring  (Apr-Jun),  summer  (Jul-Sep),  and  fall 
(Oct-Dec).  The  results  are  presented  in  terms  of  the  standard  deviation  of  the 
24  hour  temperature  change  aA^,  and  are  shown  in  Figures  48-51  corresponding 
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Figure  48.  Standard  Deviation  of  24  Hour  Temperature  Change  at  100  mb. 
(Jan-Mar  1964) 
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Figure  49.  Standard  Deviation  of  the  24  Hour  Temperature  Change  at  100  mb. 
(Apr-Jun  1964-5) 
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Figure  50.  Standard  Deviation  of  24  Hour  Temperature  Change  at  100 
(Jul-Sep  1964-5) 
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Figure  51.  Standard  Deviation  of  24  Hour  Temperature  Change  at  100  mb. 
(Oct-Dec  1064-5) 


to  the  four  seasons.  Figure  48  for  winter  shows  maximum  a  A_  values  near  and 
probably  north  of  Hawaii,  over  the  Rocky  Mountains,  and  from  flie  northeastern 
U.S.  and  eastern  Canada  eastward  across  southern  Greenland  and  south  of  Iceland 
to  Great  Britain. 

In  spring  (Figure  49)  the  values  decrease  nearly  everywhere  except  in  the  Gulf  of 
Alaska.  Maximum  values  also  occur  over  the  southern  Rockies  and  from  the 
northeastern  U.S.  across  the  North  Atlantic.  In  summer  (Figure  50),  maximum 
values  persist  in  the  Gulf  of  Alaska  and  from  the  northern  U.  S.  and  eastern 
Canada  eastward  across  the  North  Atlantic.  In  fall  (Figure  51),  peak  values  occur 
in  the  Gulf  of  Alaska  and  from  the  Great  Lakes  eastward  across  the  North  Atlantic. 

A  similar  analysis  for  southeastern  Australia  is  shown  in  Figure  52  for  the  27  day 
period  of  flight  operations  only.  It  is  interesting  to  note  that  nearly  all  the  turbu¬ 
lence  encountered  occurred  in  the  areas  of  maximum  oAp  ,  primarily  from 
Broken  Hill  and  Cobar  southeastward  to  the  coast  and  northward  along  the  coast 
through  William stown. 

For  most  areas,  some  problems  arise  in  attempting  to  establish  approximate  rela¬ 
tionships  between  aAp  and  numerical  values  for  frequency  of  occurrence  of  turbu¬ 
lence.  In  the  first  place,  the  aircraft  frequently  went  far  from  its  base  of  operations 
so  that  the  average  aAp  over  the  flight  route  might  vary  considerably  from  the 
value  at  the  operations  base.  In  addition,  flight  tracks  frequently  varied  from  one 
day  to  the  next.  This  required  estimating  an  average  value  for  aAp  for  a  relatively 
large  area  surrounding  most  bases. 

The  Australian  data  present  the  best  possibilities  of  establishing  an  approximate 
relationship  between  oAT  and  percentage  of  time  in  turbulence.  The  aAp 
values  are  fairly  uniform  over  most  of  the  flight  routes  which  correspond  closely 
to  the  areas  of  maximum  aAp  .  It  must  be  assumed,  of  course,  that  the  52  hours 
of  flight  on  11  different  days  are  representative  of  the  average  conditions  prevailing 
during  the  27  day  period  of  meteorological  data.  In  Australia,  the  observed  time 
in  turbulence  in  the  altitude  band  50-55, 000  ft  was  13%  of  the  total  flight  time  in 
this  band,  and  the  average  value  of  aAp  over  the  routes  where  the  turbulence 
occurred  is  estimated  as  3.2°C  from  Figure  52. 

The  aAp  values  in  the  Northern  Hemisphere  (Figures  48-51)  do  not  correspond 
to  the  flight  times,  but  were  derived  from  a  longer  record  so  that  a  more  nearly 
average  picture,  (but  still  not  a  true  climatological  one)  could  be  presented.  When 
the  above  relationship  between  percentage  of  time  in  turbulence  and  aAp  values 
based  on  the  Australian  data  are  extrapolated  to  the  Northern  Hemisphere,  we  con 
compare  the  observed  percentage  of  time  in  turbulence  at  each  site  with  the  longer 
term  normal.  The  results  are  shown  in  Table  V.  The  computed  values  are  based 
on  a  linear  relationship  between  aAp  and  frequency  and  using  the  Australian  data 
as  a  base,  i.e. ,  the  computed  percentages  of  time  in  turbulence  are  based  o^  the 
relationship  %  =  13  (aAp/3.2).  An  average  value  for  oAp  was  estimated  from 
Figures  48-51  for  each  operational  base  depending  upon  the  season  and  routes  of 
the  actual  flights.  In  some  cases,  the  flights  extended  for  hundreds  of  miles  from 
the  operational  base  so  that  aAT  had  to  be  estimated  for  a  fairly  large  area. 
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Figure  52.  Standard  Deviation  of  24  Hour  Temperature  Changes  at  100  mb. 
(Southeastern  Australia  17  Jul  -  12  Aug  1966) 


TABLE  V 


COMPUTED  AND  OBSERVED  PERCENTAGE  OF  TIME  IN  TURBULENCE  AND 
a £>t  VALUES  FOR  VARIOUS  SEASONS  AND  AREA 


Est.  Av.  a  At  %  of  Time  in  Turbulence 


Area 

Season 

from  Figs.  48-52 

at  50-55, 
Computed 

000  ft 
Observed 

Australia 

Winter 

3.20 

13 

Edwards,  Calif, 

Fall,  Winter 

2.50 

10 

12 

Hawaii 

Spring 

2.25 

9 

7 

Hanscom,  Mass. 

Summer,  Fall 

2.20 

9 

9 

Elmendorf,  Alaska 

Winter 

1.80 

7 

7 

Ramey,  P.  R. 

Fall 

1.50 

6 

11 

A  comparison  of  the  computed  and  observed  values  for  percentage  of  time  in  turbu¬ 
lence  in  Table  V  shows  that  they  are  in  good  agreement  for  all  areas  except  Puerto 
Rico  where  most  of  the  turbulence  was  encountered  several  hundred  miles  away, 
mostly  to  the  north  or  east.  It  is  possible  that  had  more  computational  points  been 
used  in  this  area  in  constructing  the  aA<j>  maps,  a  more  representative  value 
could  have  been  obtained.  On  the  other  hand,  it  is  also  possible  that  the  meteoro¬ 
logical  conditions  which  existed  during  the  month  of  flight  operations  were  not 
representative  of  the  longer  term  average. 

The  results  of  this  analysis  illustrated  in  Figures  48-52  and  Table  V  indicate  that 
the  standard  deviation  of  the  24  hour  temperature  change  at  HICAT  flight  levels 
may  be  a  useful  guide  for  indicating  geographical  regions  of  maximum  frequency 
of  occurrence  of  high  altitude  clear  air  turbulence.  Further  study  of  this  method 
of  approach  is  warranted. 
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SECTION  VI 


CLASSIFICATION  OF  fflCAT  BY  CHARACTERISTICS  OF  THE 
EARTH’S  SURFACE 


Monsoon  winds,  sea  breezes,  and  waves  in  the  atmosphere  on  the  lee  side  of  moun¬ 
tain  ridge  8  are  examples  of  the  association  between  macro  and  me  so  scale  meteoro¬ 
logical  conditions  and  the  characteristics  of  the  surface  of  the  earth.  In  the  regions 
of  the  atmosphere  that  are  relatively  close  to  the  earth's  surface  the  relationship 
between  the  occurrence  of  clear  air  turbulence  and  the  nature  of  the  underlying 
surface  is  relatively  well  established.  Correlations  between  HICAT  and  terrain 
features  have  also  teen  confirmed  and  the  results  are  presented  below. 


The  first  step  taken  in  the  search  for  significant  correlations  was  to  classify  the 
observations  of  turbulence  and  non-turbulence  according  to  the  characteristics  of 
the  surface  of  the  earth  in  the  region  over  which  the  WU-2  was  flying.  Desirable 
features  of  such  a  classification  system  are: 


a)  The  system  be  objective  to  the  maximum  practicable  extent 

b)  The  number  of  categories  be  small  enough  to  have  an  adequate  number 
of  samples  for  statistically  significant  results 

c)  The  boundaries  between  the  categories  be  chosen  to  produce  a  minimum 
number  of  marginal  samples. 


Crooks,  Hoblit,  and  Prophet  (2)  classified  each  turbulence  run  for  WU-2  flights 
54  to  175  according  to  a  system  which  used  four  terrain  categories.  The  assign¬ 
ment  to  a  specific  category  was  determined  by  the  character  of  the  earth's  surface 
directly  below  the  aircraft.  No  specific  cone  of  view  or  area  was  specified.  A 
detailed  study  indicated  that  more  significant  correlations  could  be  obtained  between 
turbulence  and  related  terrain  if  the  categories  were  modified  and  reduced  to  three. 
The  original  and  modified  classifications  are  as  follows,  where  h  refers  to  relief 
differences: 


Crooks,  et.  al.(2) 

Water 

Flatland  (h  <  500  ft) 
Hills  (500  <  h  <  2500  ft) 
Mountains  (h  >  2500  ft) 


Modified  System 
Water 

Flatland  (h  <  3000  ft) 
Mountains  (h  >  3000  ft) 


The  field  of  view  in  the  modified  category  was  designated  as  50  miles  upwind.  If 
a  mountain  range  was  within  50  miles  upwind  of  the  NADIR  point  of  the  aircraft 
then  the  observation  was  placed  within  the  mountain  category.  Isolated  peaks  on 
plains  or  islands  were  classified  as  flatland  or  water,  respectively.  The  application 
of  the  modified  classification  to  WU-2  flights  54  through  175  gave  results  that  are 
summarized  and  compared  to  those  of  the  initial  classification  in  Table  IV. 
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TABLE  VI 


HICAT  OBSERVATIONS  BY  TERRAIN  CATEGORY 

Number  of  Turbulent  Runs 

Terrain  Type  Crooks,  et.  al.(2)  Modified  System 


Water 

ISC 

189 

Flatland 

30 

58 

Hills 

56 

— 

Mountains 

99 

147 

A  comparison  has  been  made  betweer.  the  three  categories  of  terrain  and  the  percent 
of  turbulence  miles  per  flight  miles,  encountered  for  WU-2  flights  54  to  175.  The 
geographical  distance  covered  and  the  distance  in  turbulence  were  measured  from 
flight  track  maps  for  each  terrain  type.  Search  pattern  tracks  were  treated  as 
straight  line  distances,  thereby  eliminating  some  of  the  bias  in  data  analysis  that 
results  from  flying  several  passes  through  a  turbulent  region. 

The  influence  of  terrain  on  HICAT  is  evident  (Table  VH). 

TABLE  VH 

TOTAL  FLIGHT  MILES  (IN  HUNDREDS  OF  MILES)  FOR  EACH  TERRAIN  TYPE 
AND  PERCENT  OF  FLIGHT  MILES  IN  TURBULENCE,  FLIGHTS  54  TO  175 


Water 

Flatland 

Mountains 

Total 

>  Moderate 

Edwards 

0 

5 

(0) 

98 

(7.1) 

103 

(6.8) 

(1.2) 

Hickam 

206 

(5.9) 

0 

36 

(6.5) 

242 

(6.5) 

(0.8) 

Christchurch 

90 

(5.9) 

0 

73 

(13.0) 

163 

(9. 0) 

(2,7) 

Laverton 

1 

10) 

94 

(8.9) 

71 

(13.0) 

166 

(10.6) 

(5.4) 

Hanscom 

2 

(45.0) 

196 

(2.5) 

69 

(4.6) 

177 

(3.8) 

(0.2) 

Ramey 

132 

(5.0) 

0 

0 

132 

(5.0) 

(0.8) 

Elmendorf 

46 

(2.8) 

0 

173 

(1.9) 

219 

(2.1) 

(0.6) 

All  Bases 

477 

(5.5) 

205 

(5.4) 

520 

(6.9) 

1202 

(6.1) 

(.  16) 

>  Moderate 

(1.0) 

(1.8) 

(2.1) 

(1.6) 

The  percent  of  flight  miles  in  turbulence  is  6. 9%  over  mountains  and  approximately 
5. 5%  over  water  and  flatland.  At  five  bases  a  higher  percent  of  flight  miles  were 
turbulent  over  mountains  than  over  flatland  or  water.  Flights  from  Elmendorf, 
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Alaska,  had  more  turbulence  over  water  than  mountains  but  the  total  turbulent 
miles  were  small  (2. 1%).  On  the  other  hand,  flights  flown  out  of  Christchurch, 

New  Zealand,  and  Laverton,  Australia,  encountered  turbuL  nee  13%  of  the  distance 
over  mountains  fnd  significantly  lesser  amounts  over  other  terrain. 

The  amount  of  moderate  turbulence  over  mountains  for  all  bases  (2. 1%)  was  greater 
than  over  flats  (1. 8%)  and  water  (1. 0%).  The  figure  for  the  mountains  is  probably 
lower  than  that  which  might  be  expected  for  a  sample  covering  all  types  of  terrain 
for  all  latitudes.  In  Alaska,  where  33%  of  the  total  miles  over  mountains  were  flown, 
only  a  small  amount  of  turbulence  was  measured.  Had  some  of  those  miles  been 
flown  over  water  and  flatlands,  it  is  likely  that  the  percent  of  turbulence  would  not 
have  been  greater. 

Table  VID  presents  an  average  of  the  vertical  gust  velocity  root  mean  square 
(RMS  Uy)  values  for  each  terrain  category.  All  available  RMS  values  for  flights 
54  to  281  were  used  in  this  comparison. 

TABLE  vm 

AVERAGE  RMS  Uv  VALUES  FOR  EACH  TERRAIN  TYPE  AND  PERCENT  OF  OCCUR¬ 
RENCES  EXCEEDING  CERTAIN  LIMITS,  FLIGHTS  54  TO  281 


Water 

Flatland 

Mountains 

No.  of  Spectra 

53 

73 

107 

Avg.  RMS  Uy 

0.86 

1.02 

1.89 

RMS  Uy  >  1.00  ft/sec 

23% 

41% 

81% 

RMS  Uv  >  2.00  ft/aec 

2% 

8% 

36% 

The  effect  of  mountainous  terrain  on  turbulent  conditions  in  the  stratosphere  is 
evident,  considering  that  the  data  were  collected  over  all  four  seasons  and  from 
several  geographical  areas. 

Terrain  features  themselves  are  not  the  controlling  factors  in  high  altitude  clear  air 
turbulence.  For  example,  if  a  water  surface  were  chosen  as  one  category  then  there 
would  be  little  variation  of  clear  air  turbulence  over  the  earth  because  most  of  the 
earth's  surface  is  ocean.  In  addition,  terrain  features  such  as  high  mountains  may 
have  dissimilar  effects  upon  clear  air  turbulence  under  differing  meteorological 
conditions.  Strong  winds  perpendicular  to  the  mountain  ridge  would  probably 
yield  significantly  more  turbulence  than  light  winds  parallel  to  the  ridge.  The 
direction  of  flight  with  reference  to  a  mountain  ridge  may  have  an  important 
effect. 
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SECTION  VII 


CONCLUSIONS 


An  analysis  of  the  data  obtained  from  the  HICAT  flight  program  and  of  the  pertinent 
meteorological  data  has  led  to  the  following  conclusions: 

1.  The  power  spectra  representing  the  clear  air  turbulence  observed  over 
mountains  have  slopes  and  characteristics  that  are  distinctly  different  from 
the  slopes  and  characteristics  of  the  spectra  obtained  from  flights  over 
flatlands.  The  sample  obtained  from  flights  over  water  was  too  small 

to  justify  specific  conclusions. 

2.  No  conclusive  evidence  was  found  for  a  significant  relationship  between 
category  of  power  spectra  and  meteorological  conditions  associated  with 
HICAT. 

3.  Theoretical  analyses  indicated  definite  relationships  between  large  hori¬ 
zontal  atmospheric  temperature  gradients  and  turbulence  and  between 
relatively  large  changesln  vertical  temperature  gradients  and  turbulence. 

4»  The  turbulent  and  non-turbulent  cases  had  distinctly  different  scalar  and 
— —  vector  vertical  wind  shears  of  the  horizontal  winds. 

5.  Richardson’s  number  for  the  HICAT  levels  was  significantly  higher  than 
that  reported  in  the  literature  for  turbulent  situations  in  the  troposphere. 
However,  Richardson’s  number  did  serve  to  distinguish  between  turbu¬ 
lent  and  non-turbulent  cases  to  a  useful  extent. 

6.  Practical  forecasting  experience  for  the  HICAT  flight  program  has  indi¬ 
cated  that  strong  horizontal  temperature  gradients  and  layers  with  large 
changes  in  the  vertical  temperature  gradients  are  excellent  indicators 
of  regions  of  turbulence.  Conversely,  small  horizontal  temperature 
gradients  and  lajers  wita  uniform  vertical  temperature  gradients  indi¬ 
cates  regions  of  smooth  flying  conditions. 

7.  Relatively  large  temperature  changes  were  observed  along  the  flight  path 
when  the  aircraft  was  in  moderate  or  severe  turbulence.  During  smooth 
flight  conditions  small  or  no  temperature  changes  were  observed. 

8.  The  standard  deviation  of  the  24  hour  temperature  changes  at  the  100  mb 
level  for  the  western  part  of  the  northern  hemisphere  and  for  southeastern 
Australia  were  computed  for  each  of  the  four  seasons.  Geographical 
regions  for  which  the  standard  deviation  was  relatively  high  appear  to 
coincide  with  regions  of  relatively  high  incidence  of  turbulence.  Thus, 
this  technique  may  be  used  to  determine  the  world-wide  distribution  of 
turbulent  and  non-turbulent  regions. 

9.  Turbulence  was  observed  during  6. 9%  of  the  flight  miles  over  mountains 
and  5. 5%  of  the  flight  miles  over  water.  The  RMS  of  the  vertical  gust 
velocity  during  HICAT  was  0. 86  over  water,  1. 02  over  flatland  and  1. 89 
over  mountains. 
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